Majorana Fermions in Semiconductor Nanowires: Fundamentals, Modeling, and Experiment 



Tudor D. Stanescu^ and Sumanta Tewari^ 

^Department of Physics, West Virginia University, Morgantown, WV 26506. 
^Department of Physics and Astronomy, Clemson University, Clemson, SC 29634, USA 

After a recent series of rapid and exciting developments, the long search for the Majorana fermion - the 
elusive quantum entity at the border between particles and antiparticles - has produced the first positive ex- 
perimental results, but is not over yet. Originally proposed by E. Majorana in the context of particle physics, 
Majorana fermions have a condensed matter analog in the zero-energy bound states emerging in topological 
superconductors. A promising route to engineering topological superconductors capable of hosting Majorana 
zero modes consists of proximity coupling semiconductor thin films or nanowires with strong spin-orbit inter- 
action to conventional s-wave superconductors in the presence of an external Zeeman field. The Majorana zero 
mode is predicted to emerge above a certain critical Zeeman field as a zero-energy state localized near the order 
parameter defects, viz., vortices for thin films and wire-ends for the nanowire. These Majorana bound states are 
expected to manifest non-Abelian quantum statistics, which makes them ideal building blocks for fault-tolerant 
topological quantum computation. This review provides an update on current status of the search for Majorana 
fermions in semiconductor nanowires by focusing on the recent developments, in particular the period following 
the first reports of experimental signatures consistent with the realization of Majorana bound states in semicon- 
ductor nanowire-superconductor hybrid structures. We start with a discussion of the fundamental aspects of 
the subject, followed by considerations on the realistic modeling which is a critical bridge between theoretical 
predictions based on idealized conditions and the real world, as probed experimentally. The last part is dedicated 
to a few intriguing issues that were brought to the fore by the recent encouraging experimental advances. 
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I. INTRODUCTION 

In relativistic quantum mechanics spin- 1/2 fermions are de- 
scribed by the solutions of the Dirac equation 1 1 1, 

(z7^S^ - m)ilj = 0. (1) 

Here, 7^ (/i = 0, 1, 2, 3) is a set of 4 x 4 matrices satisfy- 
ing the anti-commutation relations {7^,7^} = 2g^^ , where 
g^^ represents the Minkowski metric and 7o7/i7o = 7^ El- 
In general, the matrices 7^ have complex elements, making 
Eq. ([T]) a set of coupled differential equations with complex 
coefficients. Thus, the general solution of Eq. ([T]) repre- 
senting the fermion field is a complex bi-spinor that is not 
an eigenstate of the charge conjugation (CC) operator. Since 
charge conjugation, i/j ^ i/j*, maps a particle into its anti- 
particle | 3 1, a complex solution of Eq. ([T]) represents a fermion 
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that has a distinct anti-fermion, with the same mass and spin 
but opposite charge and magnetic moment, as its counter- 
part. Therefore, the field of a relativistic fermion that coin- 
cides with its own antiparticle, should it exist, is necessarily 
an eigenstate of CC that must be described by a real solu- 
tion x{x) of Eq. ([T]). Real solutions of the Dirac equation are 
possible, provided one can find a suitable representation of 
the matrices 7^ characterized by purely imaginary non-zero 
matrix elements. Such a representation of the 7^ matrices 
that renders Eq. ([T]) purely real was found by E. Majorana 
in 1937 |4|. The real solutions x(^) correspond to charge- 
neutral fermions, known as Majorana fermions (MFs), rep- 
resenting particles that are their own antiparticles. Charge 
neutrality and the identification of the particle with its own 
anti-particle are not uncommon features among bosons - pho- 
tons and TT^ -mesons being two standard examples - but they 
are quite special among fermions. The neutron, for exam- 
ple, which is a charge-neutral fermion, has an anti-particle 
(the anti-neutron) that is distinguished from the neutron by 
the sign of its magnetic moment. Also, neutrinos produced in 
beta-decay are thought to be charge-neutral and have a small 
but non-zero rest mass (so that they cannot be Weyl fermions 
0, which are massless), but whether neutrinos are Dirac or 
Majorana fermions is still an unsettled question in particle 
physics. 

While Majorana fermions have remained undetected in 
high energy physics for almost 70 years, in the last decade 
they have made an emphatic entrance into the realm of con- 
densed matter physics iSHTOl. In this context, the term Ma- 
jorana fermion does not refer to an elementary particle, but 
rather to quasiparticles corresponding to collective excitations 
of an underlying complicated many-body ground state. As 
shown by N. Read and D. Green 1 11 1, in the so-called weak- 
pairing phase of a two-dimensional (2D) spinless Px + Wy 
superconductor (or superfluid), the quasiparticles satisfy an 
equation - the Bogoliubov-de Gennes (BdG) equation - that 
is similar to the Majorana form of the Dirac equation. These 
quasiparticles represent the solid state analog of the Majo- 
rana fermions from high-energy physics and, in general, are 
characterized by a finite mass and have finite energy. How- 
ever, in the context of the recent developments, the term "Ma- 
jorana fermion" has a slightly different meaning. In a su- 
perconductor, the spinless fermionic Bogoliubov quasiparti- 
cles satisfy particle-hole symmetry, 7^ = 7-^;, and, con- 
sequently, the zero-energy quasiparticles can be viewed as 
fermions that are identical with their own anti-quasiparticles, 
7o = 7o- These BdG quasiparticles with zero excitation en- 
ergy emerge as localized states bound to defects in the su- 
perconductor, where the order parameter amplitude vanishes, 
e.g., vortices and sample edges |11|. Furthermore, in two 
dimensions the quantum quasiparticles associated with these 
zero-energy bound states were shown to obey a form of quan- 
tum statistics known as non-Abelian statistics |11-15|. Ex- 
changing two particles that obey non-Abelian statistics rep- 
resents a non-commutative operation. Recently, the interest 
in the properties and the possible realization of non-Abelian 
zero-energy quasiparticles has increased dramatically, after 
they have been proposed L16J as possible building blocks for 



fault tolerant topological quantum computation (TQC) f6l. 

In this review, the terms "Majorana fermion'\ ''zero-energy 
Majorana state'\ or ''Majorana bound state'' refer to local- 
ized, charge-neutral, zero-energy states that occur at defects 
and boundaries in topological superconductors. The creation 
operator for such a zero energy state is a hermitian second 
quantized operator 7t = 7 that anti-commutes with other 
fermion operators and satisfy the relation 7^ = 1 . In a su- 
perconductor, a non-degenerate localized zero-energy eigen- 
state enjoys a form of topological protection that makes it 
immune to any local perturbation that does not close the su- 
perconducting gap. Such perturbations cannot move the state 
away from zero energy because of the particle-hole symmetry 
and the non-degeneracy condition. Since small perturbations 
to the BdG differential equation are not expected to change 
the total number of solutions, it necessarily follows that local 
perturbations (i.e., perturbations that do not couple pairs of 
MFs) leave the non-degenerate zero energy eigenvalues un- 
perturbed. This argument also implies that the zero-energy 
solutions are characterized by vanishing expectation values 
for any local physical observable, such as charge, mass, and 
spin. Otherwise, in the presence of a non-zero average, a lo- 
cal field that couples to the corresponding operator would shift 
the energy of the state. Thus, our use of the term "Majorana 
fermion" is more restrictive than that in the particle physics 
literature: the MFs in this review are charge-less, mass-less, 
spin-less, i.e., free of any internal quantum number, and obey 
non-Abelian statistics, while the MFs in particle physics sat- 
isfy standard Fermi statistics and, although charge-neutral, 
since they are eigenstates of the charge-conjugation opera- 
tor, can be massive and spin-full . Finally, as zero energy 
MFs in solid state systems are topologically protected, they 
can be removed from zero energy only by driving the system 
through a topological quantum phase transition (TQPT) ifTTl 
characterized by the closing of the energy gap at the topologi- 
cal critical point, where the MFs become entangled with other 
gapless states ifTll . 

Majorana fermions have recently been proposed in low 
temperature solid-state systems in the context of fractional 
quantum Hall (FQH) effect ETHIS [EHSl, spinless chi- 
ral p-wave superconductors/superfluids |[TT] [141 ED EH^ in 
particular in strontium ruthenate, via the realization of the 
so-called half-quantum vortices I24H261I , heterostructures of 
topological insulators (TI) and s-wave superconductors Il27l - 
I3TI , metallic surface states |32|, matallic ferromagnet-cuprate 
high-Tc superconductor 1331 , non-centro-symmetric super- 
conductors (34J |35l, helical magnets 136] , carbon nan- 
otubes |[37l , and cold fermion systems with s-wave Fesh- 
bach resonance plus artificially generated spin-orbit coupling 
and Zeeman splitting [38^, 39} . It has been also been shown 
that MFs can exist as quasiparticles localized in the topo- 
logical defects and at the boundaries of a spin-orbit (SO) 
coupled electron-doped semiconductor 2D thin film (401 ffll 
or ID nanowire ll4TH43l with proximity induced s-wave su- 
perconductivity and an externally induced Zeeman splitting. 
More recently, it was shown that similar mechanisms can 
lead to the emergence of MFs in hole-doped semiconduc- 
tor structures |i44||45l. Although the physics responsible for 
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the emergence of the zero energy MFs is identical in both 
2D and ID systems (in 3D a similar mechanism produces a 
so-called Weyl superconductor/superfluid ||46| |47l character- 
ized by gapless, topologically protected, Weyl fermion nodes 
in the bulk), the nanowire setting has some experimental ad- 
vantages, such as the option of generating the Zeeman split- 
ting by using a parallel magnetic field |4?, "43] and a signif- 
icantly enhanced gap (the so-called mini-gap) that protects 
the end-localized MFs from thermal effects |48|. In addi- 
tion, the nanowires can be arranged in networks that allow the 
maniputation of MFs for TQC |49|. The ID version of the 
semiconductor-superconductor hetero structure ifTUffi l - the 
so-called semiconductor Majorana nanowire - is a realization 
of the ID topological superconductor (TS) model first pro- 
posed by Kitaev | 23 1 in the context of TQC. In the presence of 
a small Zeeman splitting F, the Majorana nanowire with prox- 
imity induced 5-wave superconductivity is in a conventional 
superconducting state with no MFs, while for F larger than a 
critical value F^, localized MFs exist at the ends of the quan- 
tum wire. In principle, the MFs that emerge in a Majorana 
nanowire can be detected either by measuring the zero-bias 
conductance peak (ZBCP) associated with tunneling into the 
end of the wire I411I5QH541 , or by detecting the predicted char- 
acteristic fractional ac Josephson effect | 23][28j|42l|43j[55l. In 
the last year, the semiconductor Majorana wire, which is the 
focus of this review, has attracted considerable attention as a 
result of the recently reported experimental evidence for both 
the ZBCP |[56H58i and the fractional ac Josephson effect in 
the form of doubled Shapiro steps [i59i . 



II. THEORETICAL BACKGROUND 

A. Spinless {px + ipy) superconductor/superfluid 

The spinless {px + ipy) superconductor (superfluid) is the 
canonical system that supports zero energy MFs localized at 
the defects of the order parameter, such as vortices and sam- 
ple edges |11|. In 2D the mean field Hamiltonian for such a 
system is given by, 

p p 

(2) 

where, = /2m — ep with ep the Fermi energy and spin 
indices are omitted because the system is considered spinless 
(or spin-polarized). The lowest-energy solution of the BdG 
equations H^D^^ir) = E'^{r) near a vortex or near the sam- 
ple edges (where the order parameter Aq vanishes), calculated 
under appropriate conditions, is non-degenerate and has zero 
energy (E = 0), while the corresponding second-quantized 
operator (the creation operator for the Bogoliubov state) is 
hermitian, 7^ = 7. In ID the zero energy MF solutions occur 
near the ends of the wire and, as argued by Kitaev | 23 1, should 
be observable in a fractional AC Josephson effect-type exper- 
iment. Although a spinless p-wave superconductor does not 
exist in nature, this example makes it clear that the three main 
ingredients that are important for realizing zero-energy MFs 



in a condensed matter system are superconductivity, chirality 
(i.e., the + ipy orbital form of the order parameter), and 
the spinless or spin-polarized nature of the system (which en- 
sures that the zero energy solution localized at a defect of the 
order parameter is non-degenerate). 

Despite the possibilities opened by superconducting stron- 
tium ruthenate |60| and cold fermion systems in the pres- 
ence of a p-wave Feshbach resonance | 61 - 64] of realizing 
the physical conditions necessary for the emergence of MFs, 
actually realizing and observing the MFs in these systems 
are challenging tasks. In strontium ruthenate, even if the 
required half-quantum vortices can be realized, the mini- 
gap ~ A^/ei? ~ O.lmK (with ep the Fermi energy and A 
the magnitude of the p-wave superconducting order parame- 
ter) that separates the zero energy MF bound states from the 
higher energy regular BdG excitations also localized at the 
vortex cores is unrealistically small. On the other hand, in 
cold fermion systems with p-wave Feshbach resonance in the 
unitary limit, even if the mini-gap ^ A'^/ep ^ can be 
relatively large, the short lifetimes of the p-wave pairs and 
molecules represents a major experimental challenge. 

The first major attempt to reproduce the essential physics of 
spinless p-wave superconductors using only s-wave pairing 
was made in Ref. |27|, by using a heterostructure consisting 
of a 3D strong topological insulator and an s-wave supercon- 
ductor. In the cold fermion context, Ref. |38| showed that 
(also see Ref. |39|) an artificial, laser-induced, Rashba spin- 
orbit coupling and a Zeeman field, in conjunction with s-wave 
pairing interactions, can reproduce the topological superfluid- 
ity and MFs of spinless chiral p-wave superfluids using an s- 
wave Feshbach resonance. Later, the same approach that was 
used in the cold fermion context was applied to a solid state 
heterostructure where the spin-orbit coupling, Zeeman field, 
and s-wave superconductivity are provided by a semiconduc- 
tor, a magnetic field or the proximity to a magnetic insulator, 
and the proximity to an s-wave superconductor, respectively. 
The proposed setups are shown in Fig. [T] As discussed be- 
low, in the semiconductor-superconductor heterostructure the 
two main ingredients (in addition to superconductivity) that 
are essential for realizing MFs, i.e., chirality and spinlessness 
, are provided by a Rashba-type spin-orbit coupling and by 
restricting the low-energy physics of the semiconductor to a 
single (or an odd number of) relevant Fermi surfaces. 



B. Semiconductor-superconductor heterostructure 

The Rashba SO-coupled semiconductor (e.g., InSb, In As) 
with proximity induced s-wave superconductivity and Zee- 
man splitting is mathematically described by the following 
Bardeen-Cooper-Schrieffer (BCS)-type Hamiltonian: 

H = {r]k^ - ii)t^ + Fn • tr + |(A^ X tr) • ir^ + Ar^, (3) 

where r] = 1/m*, with m* being the effective mass of the 
charge-carriers, ji is the chemical potential measured from 
the bottom of a pair of spin sub-bands (in this review by 
"band", we always mean a pair of a spin-split sub-bands. 
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the band themselves being separated by the energy gaps due 
to lateral confinement), a is the Rashba spin-orbit coupling 
constant, and A is the s-wave superconducting pair-potential 
in the semiconductor, which is assumed to be proximity in- 
duced from an adjacent superconductor. In addition, n is a 
suitably chosen direction of the applied Zeeman spin splitting 
given byP = ^g/j^sB with g the effective Lande g-f actor, 
B the applied magnetic field and /j^b the Bohr magneton. 
Note that H is written in terms of the 4-component Nambu 
spinor {u^{r), ii|(r), V|(r), — and that the Pauli ma- 

trices ax^y^zj^x,y,z ^c^o^^^^ spin and particle-hole spaces, re- 
spectively. For a Zeeman splitting T larger that a critical value, 
H describes a 2D topological superconductor when the direc- 
tion of the Zeeman field is perpendicular to the plane contain- 
ing the vector k = {kxjky), i.e., for n = £. A ID topological 
superconductor with k = kx obtains for either n = ;r or 
fi = z. In general, a finite topological superconducting gap 
develops only if the component of the Zeeman field perpen- 
dicular to the effective k-dependent Rashba spin-orbit field is 
larger than a critical value. 

The Hamiltonian in Eq. ([3]) has recently been studied exten- 
sively in the context of 2D, ID and quasi- ID (multichannel) 
semiconductors Il38ti43l [65^711 . Two setups proposed for the 
experimental realization of the Majorana physics described by 
this Hamiltonian are shown in Fig. [T] The system is charac- 
terized by a topological quantum critical point (TQCP) that 
corresponds to the critical Zeeman field Tc = yA^+T^' 
where the quantity Co = (A^ + /i^ — F^) changes sign. For 
Co > 0, the (low-F) state is an ordinary, non-topological su- 
perconductor that includes perturbative effects from the Zee- 
man field and tha spin-orbit coupling but does not host MFs. 
However, for Co < 0, the (high-F) state has non-perturbative 
effects from a and can support zero-energy MF states local- 
ized at the defects of the pair-potential A. Interestingly, the 
proximity-induced pair-potential A itself remains non-zero 
and continuous across the TQCP ||4T1|72|, and, consequently, 
the two superconducting states break exactly the same sym- 
metries, namely the gauge and time-reversal symmetries. In 
the absence of topological defects and boundaries, the single 
particle spectrum of the high-F phase is similar to that of the 
low-F phase, as both are completely gapped in the bulk. How- 
ever, the high-F topological state can be distinguished from 
the non-topological superconductor at F < Fc by probing the 
topological defects and the boundaries, as they host MFs only 
for F > Fc. 

The topological quantum critical point Fc is marked by the 
vanishing of the single-particle minimum excitation gap Eo. 
By diagonalizing the Hamiltonian in Eq. ^ one obtains the 
lower-branch of the quasiparticle excitation spectrum, 

El = A^^i^^rl- 2^F2A2 + e~2r2, (4) 

where e = r]k'^ — fi and r| = F^ + a^/c^. For F near Fc, 
the minimum of is at k = 0, which corresponds to the 
minimum quasiparticle gap Eo, 

i?o = |r- VA2+/x'|. (5) 
Note that, Eo vanishes exactly at Co = 0, marking the loca- 
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FIG. 1 . (Color online) Proposed setups for realizing MFs in semi- 
conductors, (a) Two-dimensional semiconductor with Rashba SO- 
coupling proximity coupled to a superconductor and a magnetic insu- 
lator. The MFs are predicted to occur at the order parameter defects 
such as vortices. The required Zeeman field oriented perpendicular 
to the semiconductor (n = £ in Eq. ([5])) it proximity-induced by 
coupling the semiconductor to a magnetic insulator. If the semicon- 
ductor has both Rashba and Dresselhaus SO couplings, the required 
Zeeman splitting can be parallel to the surface and can be induced by 
a parallel magnetic field | 67 1. (b) One dimensional setup with Zee- 
man splitting parallel to the nanowire (n = x in Eq. ([Sj) generated 
by an external magnetic field. The end-of-wire MFs can be probed 
as zero bias conductance peaks in a local tunneling measurement, 
(c) The spectrum of the spin-split sub-bands (or the two top-most 
occupied sub-bands in a quasi- ID system) in the presence of Zee- 
man splitting. Figure adapted from Ref. L40J (panel a) and Ref. L41J 
(panels b and c). 



tion of the TQCP as a function of F, A, or /i. Furthermore, 
it has been shown that the quantity Co is the Pfaffian of the 
BdG Hamiltonian at A: = 0, (Co = Pf{H{k = 0)cr^r^)) 
[1351. The sign of Co, which determines whether the system 
has or 1 MFs, is related to the Pfaffian topological invariant 
Q (see Eq. ( pT) ) given by the product of the Pfaffians of the 
BdG Hamiltonian at k = and k = n |23|. Since the sign of 
the Pfaffian of the BdG Hamiltonian at /c = tt is generically 
positive, the sign of Co by itself determines whether the state 
of the semiconductor hosts or 1 MFs (or, in general, an even 
or odd number of Majorana bound states). The topological in- 
variant of the semiconductor-superconductor heterostructure 
is discussed in more detail in section HTdI 
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C. Topological class 

Recent work ll73lf75]i has established that the quadratic 
Hamiltonians describing gapped topological insulators and 
superconductors can be classified into 10 distinct topologi- 
cal classes that can be characterized by certain topological 
invariants. The 2D spinless chiral {px + ipy) superconduc- 
tor in the weak pairing phase (/i > 0) belongs to the topo- 
logical class D, whish is characterized by an integer Z top- 
logical invariant. This implies that, under appropriate condi- 
tions, the system can support Z number of chiral Majorana 
edge modes that remain protected against small perturbations. 
The analogous system of a 2D Rashba-coupled semiconduc- 
tor with r > Fc is also characterized by a Z invariant and sup- 
ports Z gapless chiral Majorana edge modes with dispersion 
Eedge{k) = (A//ci?)/c. Since the 2D system is in class D, does 
it automatically imply that the corresponding ID or quasi- ID 
system is also in class D? The topological class of the ID 
nanowire can be guessed from a heuristic dimensional reduc- 
tion argument. The MF mode in a ID wire (say, along the 
X-axis) can be viewed as the dimensionally-reduced version 
of the chiral Majorana edge mode of a 2D system on an edge 
parallel to the y axis. From the mathematical equivalence, 
HWdG{ky = 0) = H]P^^, where H]P^^ is the BdG Hamilto- 
nian near an end of the wire, and the property ^edge(O) = 0. it 
follows that the end of the wire supports a zero energy eigen- 
state. Since in 2D there are Z allowed chiral edge modes, 
it follows that in ID there must be, under appropriate condi- 
tions, Z zero energy MFs. Thus, the nanowire should also 
be characterized by a Z (not Z2) invariant, and, consequently, 
the purely ID Majorana nanowire should belong to the BDI 
topological class. 

To reveal the appropriate BDI classification for the 
nanowire and the associated Z invariant, it is necessary to 
identify a hidden chirality symmetry of the ID system| 48 , 76 1. 
In ID, the Hamiltonian in Eq. ([3]) anticommutes with a unitary 
operator S = Tx, where the chirality symmetry operator S can 
be written as the product of an artificial 'time reversal' opera- 
tor /C and a particle-hole transformation operator A = Tx ■ JC. 
Here, JC is the complex conjugation operator. The existence 
of all three symmetries - 'time reversal', particle-hole, and 
chirality - ensures that the Hamiltonian is in the BDI sym- 
metry class iTSlfTSl characterized by an integer topological 
invariant. Under invariance of this symmetry, the strictly ID 
nanowire can support an arbitrary integer number N {— Z) 
of protected zero energy MFs at each end. The situation is 
similar to that of a ID spinless p-wave superconductor that 
was shown fTTl to support any integer number of MFs at each 
end. However, for the semiconductor Majorana wire, the ex- 
perimentally realistic case of a quasi-YD system retains only 
an approximate chirality symmetry (the chirality symmetry is 
weakly broken by the inter-band Rashba coupling, which no 
longer commutes with S |48, 78 1) and allows multiple near 
zero energy modes at each end for values of the Zeeman cou- 
pling close to the confinement energy gap and above. The ex- 
istence of the near zero modes suppresses the gap protecting 
the MF end modes in certain parameter regimes characterized 
by r larger than the confinement energy t48| i78J . Further 



consequences of the chirality symmetry for the electrical con- 
ductance are discussed in Ref. |79|. 



D. Topological invariant 

The topological class and topological invariant of the 
semiconductor- superconductor heterostructure is analogous 
to that of a spinless Px + ipy superconductor. For the topo- 
logical invariant of the 2D spinless Px + ipy superconductor 
(which is in class D), we rewrite the Hamiltonian in Eq. Qin 
the particle-hole basis (cj^, c_k) as, 

= ^\<i^z ~^ ^x'^x^x ~ ^ykyTy^ (6) 

where we have allowed for different pair potentials A^^ , 
along the x^y directions. Writing the Hamiltonian in terms 
of the Anderson pseudo-spin vector |[8Ql (i(k) as (k) = 
(i(k).r, one can observe that in D = 2 all three components 
of d are non-zero. The topological invariant is an integer, as Z 
is the relevant homotopy group 7r2(6'^) of the mapping from 
the 2D k space to the 2-sphere of the 3-component unit vec- 
tor d = d/\d\ 111, 17|. On the other hand, in i:> = 1, the 
corresponding Hamiltonian can be made purely real (e.g., A^ 
drops out from Eq. ^ if the system is along the x-axis) and 
the vector d has only two components. Now since the /c-space 
is also one-dimensional, the topological invariant must again 
be in Z (class BDI) since 7ri(5'^) = Z. This invariant is sim- 
ply the winding number, 

^ = ^fjde{k), (7) 

where 0{k) is the angle the unit vector d makes with, say, the 
2; -axis on the x — z plane. The winding number counts the 
number of times the 2-component vector d makes a complete 
circle as k varies in the ID Brillouin zone. 

The topological invariant for the semiconductor- 
superconductor heterostructure can be defined in a similar 
way to the spinless chiral p-wave superconductor. One 
needs, however, to account for the higher dimensionality 
of the Hamiltonian matrix in Eq. ^ and the consequent 
generalization of the winding number invariant. It is clear 
that m D = 2 the Hamiltonian in Eq. ([3]) cannot be made 
real because of the complex Rashba term. In contrast, in 
D = 1 H can be made purely real, but the components of 
the (i- vector in the 4x4 Hamiltonian are themselves 2x2 
matrices. More generally, the BdG Hamiltonian of a TS 
system in D = 1, despite being real (thus preserving the 
chiral symmetry given by the operator S — Tx), can be a 
large 2N x 2N square matrix. Using the general framework 
for chiral symmetric systems Ell Ell ^ the integer topological 
invariant for this system can still be defined ^| by 
generalizing the concept of the (i- vector winding number for 
arbitrary dimensional matrices. Since, as discussed in the 
last subsection, H anticommutes with the chirality operator 
S = JC ■ A = Tx, H can be expressed as a block off-diagonal 
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matrix in a basis that diagonalizes the unitary operator 5: 



(8) 



Following Ref. |48]and|76| one can now define the variable, 

z{k) = exp(i6'(fc)) = Det{A{k))/\Det{A{k))\, (9) 
and calculate the topological invariant, 



W 



k=0 



dz{k) 



(10) 



which is an integer, including zero (W G Z). 

If the chiraHty symmetry is broken (it is weakly broken for 
quasi-lD multi-band wires with transverse Rashba coupling 
a^y 7^ 0), the number of exact zero energy MF at each end of 
the wire goes back to or 1 . The corresponding Z2 topologi- 
cal invariant is 

Q = sgn{Pf[H{k = 0)ayTy] x Pf[H{k = ^)a^r^]} , 

(11) 

and the topological class for the quasi- ID system reduces 
from BDI to D. Q = 1) represents the topologically 

trivial (non-trivial) state with (1) MF at each end of the 
wire. In the presence of the chirality invariance in ID, it 
can be shown that Q gives the parity of the integer topolog- 
ical invariant W, in analogy to the 2D case where Q for the 
semiconductor-superconductor heterostructure gives the par- 
ity 1 35 1 of the first Chern number topological invariant, which 
is an integer. 



E. Minigap 

The minigap for topological superconductors is defined as 
the energy of the first excited regular fermion state (above the 
zero energy MF state) in either the bulk or the edge of the sys- 
tem. This gap is responsible for the thermal protection of the 
MFs and, consequently, for the non-trivial physics associated 
with them, such as non-Abelian statistics, by protecting the 
Majorana bound state from mixing with regular fermion states 
above the minigap. In naturally occurring quasi-2D chiral-p 
wave superconductors, such as strontium ruthenate, even if the 
MFs are realized in the cores of the half-quantum vortices, 
the minigap S ~ A^/ei? ~ O.lmK is unrealistically small. 
In chiral p-wave superfluids, potentially realizable using cold 
fermions with p-wave Feshbach resonance, this is less of a 
problem because in this case A ~ and the minigap A^/e^? 
is of order A. In solid state hybrid structures, the minigap 
has to be large enough so that the experiments can access the 
MF physics at realistically achievable temperatures. This is 
one important experimental aspect in which the semiconduc- 
tor nanowire heterostructure has significant advantage over 
the corresponding 2D systems. In the ID system, the mini- 
gap naturally scales with the induced superconducting gap £^0 
[see Eq. Q in the nanowire, while this is generically not the 
case for the 2D systems (see, however, Ref. [.65 J). 



The reason why the minigap in the nanowires is greatly en- 
hanced over the 2D heterostructures can be understood from 
the dimensional reduction arguments. In going from the 2D 
plane to the ID wire, as the width Ly in the ^-direction is re- 
duced, the energies of the quantized chiral edge modes scale 
as l/Ly (recall that E{ky) = {A/kF)ky), i.e. as the inverse 
of the confinement size of the wire in the transverse direction. 
Since l/Ly diverges in the strict ID limit Ly 0, it implies 
that the minigap becomes arbitrarily large, in fact of the order 
of the bulk quasiparticle gap Eq. In principle, the end-state 
MFs protected by such a large minigap can be probed in local 
tunneling experiments without having to worry about other 
low energy states with energies comparable to experimental 
temperatures. 

In practice, for clean quasi- ID nanowires (L^ ^ Ly > 
Lz) with hard-wall confinement, the lowest excited state 
above the MF end state (for F > Fc, but less than the con- 
finement energy gap) has an energy that is a significant frac- 
tion of £^0 |70|. In quasi-lD wires with the Zeeman splitting 
comparable to the confinement energy gap, the minigap can 
be drastically suppressed due to the existence of multiple near 
zero energy modes at the same end, which is a manifestation 
of the approximate hidden chirality symmetry | 



E Phase diagram of semiconductor Majorana wire, sweet 
spots, and approximate chirality symmetry 

In a strictly ID wire with proximity-induced superconduc- 
tivity the Majorana-supporting topological phase emerges at 
Zeeman fields higher than the critical value |39| l4T1l 



(12) 



Consequently, for F < Fc the superconducting phase is topo- 
logically trivial, while for F > Fc the topological phase is 
realized. How is this simple picture modified in a quasi-lD 
system with multiple occupied confinement bands? If the ap- 
plied Zeeman field is much smaller than the spacings between 
the confinement bands, nothing much is changed, except that 
/i in Eq. (12) has to be understood as the chemical poten- 
tial measured from the bottom of the top-most occupied band 
(also called the Majorana band). However, the phase diagram 
is non-trivially modified when the Zeeman field is compara- 
ble with or higher than the spacings between the confinement 
bands. 

Consider a finite-width wire oriented along the x direction 
and proximity-coupled to an s-wave superconductor, with the 
interface perpendicular to the z direction. The effective spin- 
orbit field generated by Rashba coupling is oriented along the 
y direction. Two phase diagrams of this quasi- ID Majorana 
wire are shown in figures [2] and [3] The main difference be- 
tween the two diagrams is due to the different orientations of 
the Zeeman field: in Fig. [2] the field is perpendicular to the 
semiconductor- superconductor interface, while in Fig. |3]it is 
oriented along the wire. In both cases the Zeeman field is per- 
pendicular to the effective spin-orbit field generated by the 
Rashba coupling. Note that, when the magnetic field is paral- 
lel to the z axis, the normal state spectrum is characterized by 



FIG. 2. (Color online) Phase diagram as a function of chemical po- 
tential and Zeeman field oriented perpendicular to the SM-SC in- 
terface. White regions correspond to the topologically trivial phase, 
while colored regions indicate topologically non-trivial SC phases. 
The dashed lines show the location of the normal state SM sub- 
bands. Note the characteristic anti-crossings between the sub-bands 
of the transverse riy bands with opposite parity, e.g., between riy and 
ny-\-l. For a Zeeman field oriented along the wire (see Fig. [3|, these 
bands cross at the "sweet spots" (for example, the regions marked by 
red circles). Figure adapted from Ref. 171]. 

anti-crossings of the sub-bands corresponding to transverse 
Uy modes with opposite parity (see Fig. |2j red circles). By 
contrast, when the magnetic field is parallel to the wire, these 
bands cross at the "sweet spots" |69|, critical points where 
(in the absence of inter-band pairing) two topologically triv- 
ial and two topologically nontrivial phases meet. We empha- 
size that inter-band pairing stabilizes the topological phase in 
the vicinity of the sweet spots (see Fig. |3} more details about 
sweet spot physics are provided in Appendix [A|). By contrast, 
rotating the field away from the x direction favors the topo- 
logically trivial SC phase. 

The phase diagram in Fig. [3] reveals that, upon increasing 
the Zeeman field, the system undergoes a series of topological 
quantum phase transitions as it traverses successive topolog- 
ically trivial and nontrivial phases. These phases are charac- 
terized by a certain number N of low-energy modes localized 
at each end of the wire. In the thin wire limit, Ly 0, these 
modes become zero-energy Majorana bound states. In a fi- 
nite width wire, for N even, the low-energy modes acquire 
a small non-zero gap, while for N odd, TV — 1 modes be- 
come gaped and one mode (the Majorana mode) remains gap- 
less. The opening of the small low-energy gap is due to the 
inter-band Rashba coupling |48, 78], (see, for example, the 
inter-band Rashba term quyu' given by Eq. (|B4|in Appendix 
[BJ, that beaks the hidden chirality symmetry oTme purely ID 
model 148, 76i. 

Finally, we note that in the limit of weak semiconductor- 
superconductor (SM-SC) coupling, the location of the phase 
boundaries is closely related to the non-superconducting SM 
spectrum. However, the locations of the phase boundaries 
(in particular the boundary between the = and A^ = 
1 phases) depends strongly on the coupling strength ITQII . 
Specifically, from Eq. ([12]) the minimal Zeeman field nec- 
essary for reaching the topological phase in the single band 
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FIG. 3. (Color online) Phase diagram of the multiband nanowire with 
the Zeeman field oriented along the wire. Superconducting phases 
characterized by an odd (even) number N of nearly-zero energy 
modes are topologically nontrivial (trivial). The phase boundaries 
correspond to topological quantum phase transitions characterized by 
the vanishing of the bulk quasiparticle gap. The energy unit is given 
by the spin-orbit characteristic energy Ea = mef^a^ Nonuniform 
SM-SC coupling induces inter-band pairing and stabilizes the topo- 
logical phase near the sweet spots. Figure adapted from Ref. |70|. 

case is Fco = A. This critical field corresponds to the tran- 
sition between the A/" = and N = 1 phases at low val- 
ues of the chemical potential. For a multi-band system, we 
have to take into account the band dependence of the induced 
gap Ann' and the proxim ity-induced energy renormalization 
for details). Within the decoupled- 



(see section 



IIIC 



band static approximation (valid when the effective SM-SC 
coupling is much weaker than the inter-band spacing and the 
SC bulk gap), we find that the critical Zeeman field required 
for entering the topological phase with = 1 is character- 
ized by multiple minima Ten = Inn corresponding to values 
of the chemical potential at the bottom of each band n. Here, 
7nn is the effective SM-SC coupling corresponding to band 
n. Consequently, increasing the coupling strength will move 
that phase boundary toward higher Zeeman fields. In addition, 
the presence of proximity-induced inter-band pairing widens 
the topological phase near the sweet spots and further reduces 
the dependence of Fc on /i |70|. Experimentally, this may re- 
sults in an apparent insensitivity of the critical magnetic field 
associated with the emergence of the Majorana bound state on 
the chemical potential. 



G. Probing Majorana fermions I: zero bias conductance peak 
and fractional AC Josephson effect 

The simplest way to establish the possible presence of a 
Majorana bound state at the end of a spin-orbit coupled semi- 
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conductor wire with proximity- induced superconductivity is 
by performing a tunneling spectroscopy measurement. Since 
the MF state is a zero energy bound state localized at the 
end of the wire, it is expected to produce a zero bias con- 
ductance peak in the tunneling conductance dl/dV, similar to 
other zero energy boundary states in superconductors 18311841 . 
Here, / is the current from the lead to the Majorana wire and 
V is the applied voltage difference. The zero bias conductance 
peak associated with tunneling into MFs is a result of reso- 
nant local Andreev reflection of the lead electrons at the lead- 
nanowire interface |52|. In a local Andreev process, the MF 
localized at the other end of the wire plays no role, provided 
the wire is long enough so that the two MF wavefunctions do 
not overlap. If there is a significant overlap between the MF 
wavefunctions, two other processes can also contribute to the 
transport current. One process corresponds to the direct trans- 
port of electrons via the conventional (Dirac) fermion state 
formed by the two overlapping MFs (STj [85l [86l . The second 
process is the crossed Andreev reflection, in which an inci- 
dent electron from a lead coupled to one end of the wire is 
followed by a hole ejected into the lead coupled to the other 
end |87|. These two processes contribute to the non-local 
conductance or transconductance through the Majorana wire 
that we discuss in the next subsection. We note that the local 
Andreev reflection process typically gives rise to a non-zero 
sub-gap conductance even in the absence of MFs at the wire 
ends |88|. To suppress this background conductance, a gate- 
induced barrier potential may be applied at the lead-nanowire 
interface. 

The presence of a zero bias conductance peak represents 
a necessary, but not sufficient, condition for the existence 
of MFs that was theoretically proposed Ell 150^541 [TOl [89]- 
[91 1 (see Fig. |4]) and experimentally observed in recent mea- 
surements on quasi- ID semiconductor-superconductor hy- 
brid structures 156H581 . The main problem concerning the 
unambiguous identification of such a zero bias peak with end- 
localized MFs is that other sources of zero- or low-energy 
states localized at the wire-boundary f92}, "931 may produce 
similar zero-bias peaks. A more conclusive signature of MFs 
in a charge transport measurement is the quantized value of 
the T = zero bias conductance |52|. This value should 
be 2e^ /h for tunneling into a non-degnerate MF at the wire 
end, zero for the wire in the topologically trivial phase with no 
low energy end state, and 4e^ /hin the presence of a conven- 
tional (non-Majorana) zero energy state at the wire end. Such 
a conductance quantization, which can be taken as a smoking 
gun signature of MFs in semiconductor Majorana wires, was 
not yet observed and may be difficult in realistic experimental 
conditions |94j. 

In principle, a sufficient proof for the existence of MFs 
in semiconductor wires can be obtained using the Josephson 
effect. In a Josephson junction (JJ) between two topologi- 
cally trivial s-wave superconductors maintained at a phase 
difference 0, the Josephson supercurrent / is related to 
as / = Ic sin 0, where Ic is the Josephson critical cur- 
rent. In a Josephson junction between topological supercon- 
ducting wires with end-state MFs, this relation changes to 
/ /cSin^/2 I23l[28l|53|. Consequently, if the phase dif- 
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FIG. 4. (Color online) Predicted differential conductance for tunnel- 
ing into the end of a semiconductor nanowire coupled to a supercon- 
ductor. The lines represent cuts (shifted for clarity) corresponding to 
different values of the Zeeman field. As the field is increased, the 
superconducting quasiparticle gap closes at the critical Zeeman field 
Tc that marks the topological quantum phase transition. By further 
increasing the Zeeman field, the presence of a zero energy MF at the 
wire end is revealed as a pronounced zero bias conductance peak. 
In this calculation, the chemical potential is set near the middle of 
an inter-band gap, close to a sweet spot (see Fig. [Sj. For values of 
the chemical potential close to the bottom of a band, the signature 
associated with the gap closing at the TQPT is strongly suppressed, 
although the presence of a zero bias peak still reveals the MF in the 
topological phase (for details see Sec. |IVB| in particular Fig. [9]) . 
Figure adapted from Ref. L70J . 

ference (j) is controlled by a magnetic flux threaded through 
a loop containing the Josephson junction, the periodicity of 
the Josephson current changes from 27r to 47r as the semicon- 
ductor Majorana wire is tuned through the TQPT f42, 43 1. In 
practice, the doubling of the period of the Josephson effect can 
be measured by applying a small voltage V across the junc- 
tion and measuring the AC Josephson effect. The frequency 
of the AC Josephson effect in the TS phase, eV/h, is half 
the conventional Josephson frequency, 2eV/h, obtained in the 
topologically trivial superconducting phase. Note that the use 
of the AC Josephson effect is necessary to avoid the inelas- 
tic relaxation between the two branches of the local Andreev 
bound states involved in the Josephson effect li55J . 

There are several ways to understand the change in the pe- 
riodicity of the AC Josephson effect in the presence of MFs. 
Mathematically, it can be understood as a change in the de- 
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pendence of the energies of the Andreev bound states at the 
junction on the phase difference (j). While in the absence 
of the MFs the energy of the Andreev bound state goes as 
E (X coscj), with MFs this formula changes to E oc cos(j)/2. 
Since the Josephson supercurrent / is related to the Andreev 
bound state energy as / = {2e/h)dE/d(j), it follows that the 
period of the supercurrent, / ex sin(/)/2, is doubled in the 
presence of the MFs. In order for this period doubling to be 
observable, it is necessary for the two branches of the Andreev 
bound states, ±E, to maintain their occupation number as (j) 
is varied adiabatically. For DC Josephson effect this is impos- 
sible at thermodynamic equilibrium due to quasiparticle poi- 
soning effect (the higher (lower) energy state is unoccupied 
(occupied)) 1551 . In AC Josephson effect, provided the period 
of the Josephson oscillations is shorter than the inelastic re- 
laxation between the two branches ±£^(0), the Josephson pe- 
riod doubling may be observable via the fractional frequency 
eV/h. 

The doubling of the Josephson period in (/) can also be un- 
derstood from purely topological considerations. The fermion 
parity of the ground state of the topological superconductor 
(which is in the form of a ring with a weak link in the Joseph- 
son set up) flips as the threaded flux changes by 27r. With a 27r 
flux change the superconductor is thus in an excited state with 
an electron from the ground state ejected to the junction and 
held between the two MF states. Only with a total flux change 
of 47r the ground state Fermion parity flips back and the elec- 
tron can leave the state formed by the MFs. This makes the 
junction return to the original state and consequently doubles 
the period of the Josephson effect from 27r to 47r 1 1111231. 

Even though the fractional Josephson effect is a robust sig- 
nature of MFs, such an effect cannot be ruled out in ballistic 
S-N-S junctions and in JJs made of ID p-wave supercon- 
ducting wires such as the quasi-lD organic superconductors 
1551 . However, in fermion parity protected superconductors, 
the MF mediated fractional Josephson effect is topologically 
protected while there is no such robustness (for instance, to 
disorder in the ballistic S-N-S junctions, or to magnetic fields 
in p-wave wires) in the other two systems. Nevertheless, the 
fractional Josephson effect in semiconductor Majorana wires 
is susceptible to the quasiparticle poisoning effect [551, as 
well as to the non-adiabaticity effects 1951 . The latter effect, 
in particular, shows that the fractional Josephson effect cannot 
be taken as an unambiguous signature of MFs because half or 
other fractional frequencies are in principle possible even in 
junctions between conventional superconductors. Recently, 
experimental evidence of fractional ac Josephson effect in the 
form of doubled Shapiro steps has been reported in semicon- 
ductor Majorana wires |59|. In this review we will not further 
discuss this effect, but will focus instead on experiments tar- 
geting the zero bias conductance peak. 



H. Probing Majorana fermions II: quantum non-locality, 
transconductance, and interference 

Neither the non-quantized zero bias conductance peak, nor 
the fractional AC Josephson effect could constitute a suf- 



ficient proof for the existence of MFs in semiconductor- 
superconductor heterostructures, as in principle they can arise 
even in the absence of MFs. However, a conclusive experi- 
mental proof could be obtained by making use of the intrin- 
sic non-local properties of these states. In MFs, non-locality 
stems from the absence of an occupation number associated 
with them individually. To define the electron occupation 
number in terms of spatially separated MFs, one must con- 
sider a pair of MFs, 7a, 76, and define the second quantized 
electron creation operator as (i^ = 7a + ijb- The quantum 
state of the system is then determined by the eigenvalues of 
the electron occupation number operator Ud = d)d = 0, 1. 
Since Ud is related to the MF operators by 

rid ^ , (13) 

it follows that the state of the whole system is determined by 
non-local correlations between the spatially separated MFs 7a 
and 75. An idea to probe this non-locality involves injecting 
an electron into one end of the Majorana wire and retriev- 
ing it at the opposite end. By connecting leads to the left 
and the right ends of a wire that hosts MF states localized 
near the two ends, one could imagine that an electron injected 
into the end a flips the occupation number Ud from Ud = 
to Ud = 1. The injected electron can then escape from the 
end 6, flipping the occupation number of the nanowire state 
back to = 0. Such a process, where an electron can enter 
from one end and exit at the other end as an electron, can be 
viewed as Majorana-assisted electron transfer. As has been 
discussed before ||5T][85l|86l, such a transfer should not vio- 
late the locality and causality principles. In other words, in the 
absence of an overlap between the two MF wavef unctions, the 
probability of an electron appearing at end h becomes com- 
pletely independent of the lead-wire voltage difference ap- 
plied at the opposite and |51|. Hence, there is no violation 
of causality. The non-local conductance (transconductance), 
given by dli^/dVa, can only be non-zero when there is a fi- 
nite overlap between the two MF wavef unctions. However, in 
addition to the Majorana-assisted electron transfer, a second 
source 1 87 1 for non-zero transconductance is represented by 
the crossed Andreev reflection (CAR), also known as Cooper 
pair splitting. In this process an electron injection into a Ma- 
jorana bound state at one end is followed by the emission of 
a hole from a second Majorana state at the other end. The 
net result is the injection of a Cooper pair in the topological 
superconductor nanowire. The CAR contribution to transcon- 
ductance, which is opposite in sign to the contribution aris- 
ing from Majorana assisted electron transfer, is also non-zero 
only when the wavefunctions of the two MF states overlap 
ISTl . It has been shown that I961, in the case of symmetric 
tunneling between the two leads and the MFs (ta = h = t), 
the transconductance in the semiconductor Majorana wire is 
given by, 

dh^. 32K 

dVa 16T^ + {S^-V^y+8T^{5^+V^)' ^ ' 

Here, Va is the voltage at end a, lb is the current at end b, 
r oc is the lead-induced broadening of the MF level, and 
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5 is the overlap integral between the two MF wavef unctions. 
Note that the transconductance vanishes in the Hmit 6^0, 
which is consistent with earher results (51] |86j |87l. Even 
when the direct wavefunction overlap of the MFs is vanish- 
ingly small, an effective coupling between the MFs, and con- 
sequently a non-zero Majorana assisted electron transfer am- 
plitude, may be present if the topological superconductor has 
an appreciable charging energy [97 1. In this case, both the 
Majorana assisted electron transfer and CAR can result in a 
non-zero transconductance between the leads. 

Even the MF induced transconductance, while interesting 
and nontrivial, cannot be considered as a definitive signature 
of MF modes because conventional near-zero energy states 
(such as those produced by localized impurities) trapped near 
the contacts with the leads can also produce such non-local 
signatures in the presence of superconductivity [96] . On the 
other hand, a non-local tunneling spectroscopy interference 
experiment, similar to earlier interference based proposals in 
topological insulators and superconductors |29, 30, 97, 98], 
has recently been suggested as capable to provide a direct 
verification of the non-local physics of ID wires containing 
end MFs [[96|. The effect requires non-local fermion par- 
ity 1 99 1, which is unique to topological systems and cannot be 
emulated by conventional near-zero-energy Andreev states or 
any other local excitations near the wire ends. The proposed 
scheme, where the tunneling amplitude is measured as a flux 
dependence of an energy level and the fermion-parity is fixed 
by a superconducting single-electron transistor configuration, 
is suitable for unambiguous experimental testing of the pres- 
ence of MFs in semiconductor Majorana wires. 



potentials. To account for the complex phenomenology of a 
Majorana-supporting structure and make connection with ex- 
perimental observations, it is key to incorporate these factors 
in the theoretical model. The general form of a Hamiltonian 
that describes a semiconductor-superconductor (SM-SC) hy- 
brid structure is 

^tot = ^SM + ^int + i^Z + + ^SC + ^SM-SC, (15) 

where JYsm is a non-interacting model for the semiconductor 
component (or, in general, other spin-orbit coupled material, 
e.g., topological insulator), i^int contains many-body elec- 
tronic interactions, Hz describes the applied Zeeman field, 
Hy contains terms that account for disorder and gate poten- 
tials, i^sc is the Hamiltonian for the superconductor, and 
^SM-sc describes the semiconductor-superconductor cou- 
pling. Here, we will not address the problem of electron- 
electron interaction and will assume that i^int = 0. However, 
we emphasize that Coulomb interaction plays an important 
role in low-dimensional systems, as it generates a density- 
dependent renormalization of the electrostatic potential. Con- 
sequently, the Zeeman field F and the chemical potential /i 
of a Majorana wire are not independent variables 1 1001 and, 
for a given configuration of external gate potentials, the sys- 
tem is characterized by a specific function ji = /i(r). This 
has important experimentally observable consequences 1 1 001 , 
as we will discuss in Section [Tvl At the Hartree level, the 
effects of Coulomb interaction can be incorporated using a 
self-consistent scheme similar to that described in Ref. IIOII. 



Below, we derive the low-energy effective model for the het- 
erostructure described by Eq. 
i^int = 0. 



15) under the assumption 



III. REALISTIC MODELING OF SEMICONDUCTOR 
MAJORANA WIRES 



A. Tight-binding models for semiconductor nanowires 



The construction of the effective low-energy model for a 
Majorana hybrid structure involves three main steps: i) devel- 
oping a tight-binding model for the component that provides 
spin-orbit coupling (e.g., semiconductor, topological insula- 
tor, etc.) and projecting onto a reduced low-energy subspace, 
ii) incorporating the superconducting proximity effect, and iii) 
defining an effective Hamiltonian based on a linear approxi- 
mation for the frequency-dependent proximity-induced self- 
energy. Below, we provide the relevant details and point out 
the main approximations involved in this construction. In ad- 
dition, we discus several aspects of Majorana-supporting hy- 
brid nanostructures that are not captured by simple models 
of ideal ID systems, but represent critical components of the 
experimental realizations of Majorana nanowires, such as the 
presence of disorder and smooth confining potentials. 

The existence of Majorana fermions in an ideal topologi- 
cal SC system does not depend on the details of the Hamilto- 
nian. However, under realistic laboratory conditions, the sta- 
bility of the Majorana mode and the low-energy phenomenol- 
ogy of the heterostructure depend critically on a large num- 
ber of specific parameters, including details of the electronic 
band structure, nanostructure size and geometry, coupling at 
the interface, type and strength of disorder, and applied gate 



The first element in the development of a low-energy ef- 
fective theory for a Majorana-supporting hybrid structure is a 
tight-binding model for the spin-orbit coupled component of 
the system (e.g., electron-doped SM, hole-doped SM, topo- 



logical insulator, etc.) corresponding to H^m in Eq. ( p3] ). In 
general, when choosing the model, one has to strike a balance 
between accuracy and simplicity and one has to consider two 
key aspects: i) the dimensional reduction from three dimen- 
sions (3D) to quasi-2D or quasi-lD, and ii) the projection 
onto a reduced low-energy subspace. To illustrate some of 
the possible issues, we consider the case of electron-doped 
and hole-doped semiconductors. 

The simplest tight-binding model for an electron-doped 
wire (or a thin film) is a two-band model with nearest neigh- 
bor hopping. The natural starting point for constructing such 
a model is a 3D Hamiltonian for conduction electrons in the 
effective mass approximation. 



%(fc) = ^ ( 



2m* 



(16) 



where c^cr is the annihilation operator for a particle with wave 
vector k and spin a, m* the effective mass of the conduction 
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band, and /i is the chemical potential. In a thin film, the mo- 
tion in the transverse direction is quantized and the transverse 
modes can be obtained by solving the quantum problem in- 
volving the single-particle Hamiltonian corresponding to Eq. 
( p^ with the substitution kz —idz and the confining po- 
tential V{z). In the presence of a transverse field that breaks 
inversion symmetry, the spin and orbital degrees of freedom 
are coupled and the effective Rashba-type spin-orbit interac- 
tion (SOI) is described by 

'Hsoi(^) = ^R^cl{kyax - kxay)ck, (17) 



where k = (k^^ky), aR is the Rashpa coefficient, are Pauli 
matrices, we have used the spinor notation = (c^^, c^^), 
and we have assumed that only one transverse mode is rele- 
vant to the low-energy physics. The physics described in the 
long wavelength limit by the Hamiltonian Hq + Hsoi + 
where Hy represents the contribution from the confining po- 
tential, can be also determined using a tight-binding model 
defined on a lattice. For example, considering a simple cubic 
lattice with lattice constant a and nearest-neighbor hopping, 
the semiconductor Hamiltonian reads 



Hsu = i^o + Hsoi = -to 4^a^j^ - ^^Yl ^1^^^^ 

i,S,a i,a 
■ ~ + h.c] , (18) 



y 



where Hq includes the first two terms and describes nearest- 
neighbor hopping on a simple cubic lattice and the last term 
represents the Rashba spin-orbit interaction. Here, i = 
i^xj^^^z) labels the lattice sites, S G {Sx^Sy^Sz} are 
nearest-neighbor position vectors. The confining potential in- 
volves the additional local term Hy = ^ V{i)c\^Cia. The 
parameters of the tight-binding Hamiltonian ([18]) are deter- 
mined by the condition that i^sM and Ho + 'Hsoi describe 
the same physics in the long-wavelength limit k ^ ^. Con- 
sequently, the hopping parameter in ( 18 ) is to = fi^a~^/2m* 
and the Rashba coupling is a = a^^/a. 

For hole-doped SMs, an effective tight-binding Hamil- 
tonian can be obtained using a similar approach and start- 
ing from the Luttinger 4-band model ifTMl fTOSI , Ho(^) = 
c^[/il(^) — lA^k, where Ch and cj[. are four-component 
spinors. The single-particle Luttinger Hamiltonian is 



hL{k) = 



mo 



271 + 572 



(19) 



72 Y ^^'^i ~^^Y1 ^i^jJiJj 



where mo is the free electron mass, k — (^x? ^2/' ^2;) the 
wave-vector, J = ( J^;, J^, J^) is a set of 4 x 4 matrices rep- 
resenting angular momentum 3/2, and 71, 72, and 73 are the 
Luttinger parameters. We assume that the dominant contribu- 
tion to the valence band comes from p-like orbitals, 
and and that in the presence of SOI the top valence band 




FIG. 5. (Color online) Comparison of 3D spectra corresponding to 
the 8-band Kane-type model (blue lines), the 4-band Luttinger-type 
model given by Eq. ( [2T| (orange), and the 2-band model described 
by Hq in Eq. ([18]) (yellow). Each band is double spin-degenerate. 
Note that the top valence band (heavy hole band) has the same disper- 
sion in the 8-band and 4-bad models, while quantitative agreement 
between the light-hole bands and between the conduction bands is 
present only in the vicinity of the F point, i.e., for ^ 0. The 
numerical values of model parameters correspond to InSb. 

corresponds to the eigenstates of the total angular momentum 
momentum J with j = 3/2, 



-^(T|x)-z|y))x^ + y^|z)x±, 



1 

71 



(T|x)-z|r))x±, 



(20) 



where m = ±1/2, ±3/2 are the eigenvalues of Jz and x± are 
spin eigenstates with Sz = ±1/2. The effective tight-binding 
model is constructed using the eigenstates ( [2Q| ) as a basis. 
More specifically, we consider an fee lattice and nearest- 
neighbor hopping between the states <l>m(^) and 
where i and j are nearest neighbor sites, with hopping matrix 
element . The tight-binding Hamiltonian has the form 



mm' t 
ij in 



^E4 



(21) 



where ct^ is the creation operator for the state <^>rn(^)- If we 
considering the xy plane perpendicular to the (0, 0, 1) crys- 
tal axis, each site has four in-plane and eight out-of-plane 
nearest-neighbors. The corresponding hopping matrix ele- 
ments can be expressed in terms of three independent parame- 
ters, ^1,^2, and ts. For example, the in-plane diagonal compo- 



nents are t 



±3/2±3/2 



= ti/2±2t2andt;^^/^^^/^ =ti/2-2t2, 
while the diagonal out-of-plane hoppings can be written as 
^±3/2±3/2 ^ ^^/2 - t2, and Cr^'/^ = h/2 + t2. Finally, 
the values of the independent parameters are determined by 
the condition that the low-energy, long-wavelength spectrum 



of the lattice Hamiltonian pT) be identical with the spectrum 
of the Luttinger model. Explicitly, we have ti = ^yiTi^ /moo? , 
i = 1,2,3, where a is the lattice constant. For a hole-doped 
thin film with structural inversion asymmetry, the Rashba- 
type spin-orbit coupling is modeled phenomenologically by 
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adding a term similar to i^soi in Eq. (18), but for spin 3/2, 



i.e., with (T ^ J. This SOI term induces a spHtting of the 
top valence band (the heavy hole band) that is proportional to 
in the limit /c ^ 0. This is contrast with the linear split- 
ting that characterizes the Rashba splitting of the conduction 
band modeled by the Hamiltonian i^sM from Eq. ([18]). Note 



that for both the 2-band model ( 18 ) and the 4-band Luttinger- 
type model, the effective Rashba coefficient a is an indepen- 
dent parameter that does not depend on the confining potential 
V{i). Typically, in the numerical calculations V{i) is a hard- 
wall potential ( zero inside a finite region and infinite other- 
wise) that does not beak the structural inversion symmetry. 



Expanding the basis (20) to include the eigenstates of the 
total angular momentum momentum J with j = 1/2 (corre- 
sponding to the split-off band), as well as s-type states (cor- 
responding to the conduction band), allows to construct an 8- 
band tight-binding model for both electron-doped and hole- 
doped SMs that is equivalent in the long-wavelength limit 
to the 8-band Kane model 1 1041 . There are two key differ- 
ences between this model and the simpler models described 
above. First, the Rashba-type spin-orbit coupling in the 8- 
band model is determined by implicitly by the asymmetric 
confining potential V{i) and does not involve any additional 
independent parameter. Consequently, this model captures 
the correlation between the strength of the spin-orbit cou- 
pling and the transverse profile of the wave function. In turn, 
this profile plays a key role in the SC proximity effect, as we 
will show in section |IIIB[ By contrast, these correlations are 
not captured by the simplified models. The second key dif- 
ference becomes manifest when describing low-dimensional 
nanostructures. In the large length scale limit, the 2-band and 
4-band models described above generate conduction and va- 
lence band spectra, respectively, that agree quantitatively with 
the results obtained using the 8-band model. This is illus- 
trated in Fig. [5] which shows a comparison between the 3D 
spectra obtained using these models in the absence of spin- 
orbit coupling. The large discrepancies at wave vectors away 
from the F point suggest that for systems with reduced di- 
mensionality theses models will generate significantly differ- 
ent results. This behavior is illustrated in Fig. [6j which shows 
a comparison between the conduction band spectra of a SM 
film of thickness Lz = 50nm obtained using the 8-band and 
the 2-band models. In essence, the spectrum of the 8-band 
model is characterized by strong non-parabolicity and by a 
sub-band-dependent effective mass. Also, the quasi-2D ef- 
fective mass predicted by the 8-band model is strongly de- 
pendent on the film thickness. These features cannot be re- 
produced by the simplified 2-band model. Moreover, similar 
discrepancies characterize the valence bands obtained using 
the 8-band and 4-band models. All these differences become 
even more pregnant in quasi-lD nanostructures. Nonetheless, 
since Majorana physics is mainly controlled by a relatively 
small number of low-energy states, the simplified models can 
provide a reasonably accurate description of the semiconduc- 
tor system with a proper choice of effective model parame- 
ters and within a limited range for the control parameters, i.e., 
chemical potential and Zeeman splitting. 

An important aspect that has to be addressed when solving 
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FIG. 6. (Color online) Comparison between the conduction band 
spectra of a 50nm InSb film obtained using the 8-band (blue) and 
the 2-band (yellow) models. Note that both models predict the 
same value of the effective mass for the lowest energy confinement- 
induced sub-band, but the effective masses of higher energy sub- 
bands, as well as the values of the energy gaps between sub-bands, 
are significantly different. 

numerically a specific model for the Majorana wire is rep- 
resented by the large number of degrees of freedom in the 
problem. We emphasize that Majorana fermions are zero- 
energy bound states localized near the ends of a SM wire that 
is proximity-coupled to a superconductor, hence we are inter- 
ested in modeling a finite quasi- ID system with no particular 
symmetries, in the presence of disorder and external poten- 
tials. For a typical electron-doped SM nanowires modeled 



using the 2-band model ( 18 ) the number of degrees of free- 



dom is of the order 10^-10^. One possible solution is to use 
a coarse-grained lattice model with an effective lattice con- 
stant tteff much larger the actual SM lattice constant. For ex- 
ample, choosing aeff = 40a reduces the number of degrees 
of freedom to about 10^-10"^. However, while this is numer- 
ically convenient, it becomes difficult to address the short- 
range properties of the system, for example the effects of cer- 
tain types of disorder. More importantly, the properties as- 
sociated with Majorana physics are basically controlled by a 
reduced number of low-energy degrees of freedom, hence it 
is more natural to project the Hamiltonian onto the relevant 
low-energy sub-space than to perform an overall energy- 
independent reduction of the Hilbert space. In addition, the 
simple models used in the calculations are expected to be 
highly inaccurate at high energy. The technical details associ- 
ated with the projection onto a low-energy subspace are pre- 
sented in Appendix|B]for the case of the 2-band tight-binding 
model. 



B. Proximity effect in semiconductor 
nanowire-superconductor hybrid structures 

A critical ingredient of any recipe for realizing Majorana 
fermions in a solid state hybrid structure is the proximity- 
induced superconductivity. In essence, the electrons from the 
SM nanowire acquire SC correlations by partly penetrating in 
the nearby bulk s-wave superconductor. In conjunction with 
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spin-orbit coupling, these correlations generate effective p- 
type induced superconductivity in the SM nanowire. More 
specifically, each double degenerate confinement-induced 
band generates a combination of Px + ipy and Px — ipy su- 
perconductivity. Further, the Zeeman splitting breaks time- 
reversal symmetry and selects one of these two combinations, 
making the SM nanowire a direct physical realization of Ki- 
taev's toy model for a ID spinless p-wave superconductor 
||23]| . While induced pairing is the most prominent aspect of 
the proximity to the bulk superconductor, another effect is the 
renormalization of the energy scale for the nanowire. Quali- 
tatively, this renormalization can be understood as resulting 
from the reduction of the quasiparticle weight of the low- 
energy SM states due to the partial penetration of the cor- 
responding wave functions into the SC. As a consequence, 
the low-energy effective model for the nanowire will con- 
tain rescaled values of the original parameters (e.g., hopping 
parameters, spin-orbit couplings, Zeeman field, etc.), as we 
show below in section JlTCj 

To address quantitative aspects of the SC proximity ef- 
fect, one has to consider specific models for the relevant 
terms in the total Hamiltonian ( p3] ): Hsc, which describes s- 
wave bulk superconductor, and i^sM-sc, representing the SM 
nano wire-superconductor coupling. At this point, one possi- 
ble approach is to treat the pairing problem self-consistently 
I105H107 1 in order to account for spatial variations of the or- 
der parameter near the interface and effects due to interactions 
inside the nanowire IIQ8HIIIII and the presence of a mag- 
netic field. In this work, we will not discuss this aspect of 
the proximity effect and will model the bulk semiconductor at 
the mean-field level using a simple tight-binding Hamiltonian 
characterized by a constant pairing amplitude Aq. Explicitly, 
we have 

Hsc = V {ttj - fiscSij) at^a^v+Ao V(4t4;+«*i«n). 



where i and j label SC lattice sites, is the creation operator 
corresponding to a single-particle state with spin a localized 
near site i, and fisc is the chemical potential. Similarly, we 
can write the SM-SC coupling as 



SM-SC 



*0,jo ^)Cr 



h.c. 



(23) 



where io = (ix^h^ioz) and jo = {jxjyjoz) label lattice 
sites near the interface in the SM and SC regions, respec- 
tively, m is the quantum number that labels the SM states 
(e.g., if the SM is described by a 2-band model m = a), and 
i^j^ are coupling matrix elements between SM and SC local 
states. Since we are interested in the low-energy physics of 
the SM nanowire, it is convenient to integrate out the SC de- 
grees of freedom and define an effective action for the wire. 
In the Green function formalism, this amounts to including 
of a surface self-energy contribution in the SM Green func- 
tion l4nr7nfTT2ll . The basic structure of this self-energy con- 
tribution is illustrated by the local term 



Aoayi 



<rz 



(24) 



where |t| represents a measure of the SM-SC coupling, vp 
is the surface density of states of the SC metal at the Fermi 
energy, cf\ and t\ are Pauli matrices associated with the 
spin and Nambu spaces, respectively, and C is a proximity- 
induced shift of the chemical potential. Note that [24] con- 
tains an anomalous term proportional to Ty that describes the 
proximity-induced pairing. Also, the diagonal term that is 
linear in frequency in the limit ^ is responsible for the 
reduced quasiparticle weight and the corresponding energy 
renormalization. Since the basic aspects have been addressed 
by several authors ||4T] [65] [71] I112H115I , we focus here on 
two problems that are critical to understanding the proximity 
effect in real SM-SC hybrid structures: i) the role of multi- 
band physics, and ii) the dependence on specific features of 
the SM-SC coupling matrix elements. 

Proximity effect in multiband nanowires 

The key feature that differentiates the multiband case from 
its single band counterpart is the emergence of proximity- 
induced inter-band pairing. In addition, the proximity- 
induced renormalization of the SM energy scales has a ma- 
trix structure, rather than being described by an overall factor. 
To illustrate these features, we consider an electron-doped 
SM nanowire with rectangular cross section and dimensions 
Lx ^ Ly ^ Lz in contact with an s-wave SC. The nanowire 
is modeled using the 2-band mode l (p^ w ritten in the basis 
Verier (^), as described by equations (Bl -B6), and the SM-SC 
coupling is given by Eq. ( [23] ). For simplicity, we assume lat- 
tice matching between the SM and the SC and nearest neigh- 
bor hopping across the interface, t^J^ = i{iy)Si^ j^, where 
d = (0, 0, is a nearest-neighbor position vector. The de- 
pendence of i = i{iy) reflects the possibility of nonuniform 
SM-SC coupling across the nanowire, as in the recent exper- 
iments by Mourik et al. |56|. After integrating out the SC 
degrees of freedom, the effective surface self-energy can be 
(22) written in the spinor basis ipri (V^rit^ V^n,;, V^nt^ i^niV as 



E 



(25) 

where jo = io -\- d and Jq = Iq -\- d and Gsc is a matrix 
that contains both normal and anomalous terns that represents 
the Green function of the superconductor at the interface. We 
note that there are three sources of proximity-induced inter- 



band coupling in Eq. ( [25] ): i) the position dependence of Gsc, 
ii) the non-vanishing coupling between states with arbitrary 
values of and n'^ , and iii) the position dependence of i. 

The first source has not yet been explored, but could have 
significant effects when the SC itself has small character- 
istic length scales. To account for these effects, the SC 
Green function has to be calculated explicitly by taking into 
account all relevant details, including the size and geome- 
try of the system. On the other hand, assuming a large 
superconductor with a planar surface, the SC Green func- 
tion becomes Gsc = Gsci^^'^x — '^x^^y ~ '^'y) 
be expressed (TOj [71] 11121 in terms of its Fourier transform 
Gsc{^,k\\) ^ -i/F[(^+Aocr^r^)/v'Ag -uj^^Cr^]. In this 
case, since the dependence of ^5-0(^,^11) on the in-plane 
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wave-vector is very weak f9l [TOl fTT2ll . the surface self- 
energy contribution is practically local and does not represent 
an additional source of inter-band coupling. 

The second mechanism, which couples different 
confinement-induced bands, is due to the presence 
of the interface and and does require any in-plane spatial in- 
homogeneity. Qualitatively, this coupling can be understood 
in terms of virtual processes in which a electron occupying 
a state from the band tunnels into the superconductor, 
then returns to the SM wire into a state from a different 
band, n^. Assuming uniform hopping i across the interface 
, the effective coupling due to such processes depends on 
the values of the wave function at the interface and can 
be expressed as 7n,n' = ^fI^P^u, (^0;^)^n; (^o^), where 
is given by Eq. (BT|. We note that the effect of this 



inter-band coupling becomes negligible in the limit of strong 
confinement, when 7n^n' is much smaller than the inter-band 
gap. 

The third source of proximity-induced inter-band coupling 
is due to non-homogeneous interface hopping, i = i{iy). This 
mechanism is discussed in detail in Ref. |70|. The effective 
SM-SC coupling of a SM nanowire with rectangular cross 
section and nonuniform interface tunneling can be written as 

Ny 

Inn' = Sn^n'^T^F<i)nA^0z)(t)n'^{i0z) ^ l^^n^ (^2y)0n; (^2/ 

iy=l 

(26) 

Neglecting the spatial dependence of the SC Green function 
at the interface, the proximity-induced effective self-energy 
becomes 



(27) 



where 7^71,' is given by Eq. ( 26 ) and C is a constant that de- 
pends on the details of the SC band structure. Note that Eq. 
(27) is valid for frequencies inside the SC gap, \uj\ < Aq. 



Also, we emphasize that the discussion leading to Eq. ( 27 ) 
was based on the 2-band tight-binding model for the SM 
wire, which is constructed s-type localized orbitals using as 
a basis, with spin being the only internal degree of free- 
dom. In more complex models, the local states are labeled 
by a quantum number m a, or by a set of quantum num- 
bers. The hopping matrix elements across the interface, which 
parametrize the coupling Hamiltonian ( [23] ), depend explic- 
itly on these quantum numbers and, in turn, the proximity- 
induced surface self-energy will depend on the details of this 
coupling. An example illustrating the the proximity effect for 
hole-doped SM nano wires is given in Appendix [C] 



C. Effective low-energy Bogoliubov-de Gennes Hamiltonian 
for the Major ana wire 

In the presence of proximity-induced superconductivity, 
the low-energy physics of the SM nanowire is described by 
the Green function matrix Gnn'i'^) that includes nonzero 



anomalous terms. Using the results of sections III A and III B 
one can write the inverse of the Green function matrix as 



[G ^]nn'{^) =UJ - Hn 



(28) 



where i^nn' is the effective low-energy Hamiltonian for 
the SM nanowire in the Nambu space and T^rm'i^) is the 
proximity-induced self-energy. For an electron-doped wire 
described using the 2-band tight-binding model, the self- 
energy is given by Eq. ([27]), while the low-energy Hamil- 



tonian can be obtained by expanding Eq. ( |B6| ) to include both 
particle and hole sectors. Explicitly, we have 



Hjin' \^n H~ TcTq;;] T^Sfq^ji' 
-\- iddn^n'^ Qnxn'^^y^riyn'y Qnyn'y^x^z^nxn'^ 



, (29) 



where the relevant quantities are the same as in Eq. ( |B6| ), 
is a Pauli matrix in the Nambu space, and the identity ma- 
trices in the spin and Nambu spaces, ctq and tq, respectively, 
have been omitted for simplicity. The eigenvalues of the low- 
energy states can be obtained by solving the Bogoliubov-de 
Gennes (BdG) equation 



det[G-\uj)] =0, 



(30) 



where the frequency is restricted to values inside the bulk SC 
gap, I a; I < Aq, and the Green function is given by equations 
([2S]), ([ST]), and ([29]). 

Solving Eq. ( 30 ) numerically can be rather demanding, but 



one can further simplify the problem by noticing that the rele- 
vant energy scale for majorana physics (e.g., the induced pair 
potential A) is typically much smaller than the bulk SC gap 
Aq. Consequently, one can focus on the the low-frequency 
limit \uj\ <C Aq and consider the self-energy within the static 

. In the static approxima- 



approximation \/aJ 



tion, Eq. (|28l) 

"Jnn' {^y'Ty 



becomes [G 
C/AqTz), where Q 

r 



UjQr 



— ^nn' ~r ^nn' /Ao ac- 
counts for the proximity-induced normalization of the SM en- 
ergy scales. Since Qnn' is a positive definite matrix, the BdG 
equation ([30|) can be rewritten in the static approximation as 
det[cj — H^] = 0, where H^^ is a frequency-independent 
quantity that can be viewed as an effective BdG Hamiltonian 
for the SM nanowire with proximity-induced superconductiv- 
ity. Explicitly, the effective Hamiltonian has the form 



eff 



^l^nn'Tz' (31) 



The first term in the right hand side of Eq. ( [3T| ) represents 
the SM Hamiltonian renormalized by the proximity effect. 
Note that, in general, this proximity-induced renormalization 
is described by a matrix, rather than an overall factor Z repre- 
senting the reduced quasiparticle weight. The renormalization 
matrix is the solution of the equation 



(32) 



where ^mm' is the effective SM-SC coupling matrix. In the 
single-band limit, we have = Z = (1 + 7/Ao)~^. Note 
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that the renormaUzed Hamiltonian is written in the symmetric 
form ZRZ, rather than ZR, to ensure its hermiticity in the 
general, multi-band case. The second term Eq. ( [3T] ) contains 
the proximity-induced pairing matrix 



(33) 



In the single-band limit, this reduces to the induced SC pair- 
potential A = 7A0 /(7+ Aq). Finally, the last term in ( 31 ) de- 



scribes proximity-induced energy shifts and inter-band cou- 
plings that have only a rather limited quantitative relevance. 

We emphasize here a key aspect of the proximity ef- 
fect in multi-band SM-SC hybrid structure: the proximity- 
induced inter-band coupling and the induced inter-band pair- 
ing. These effects depend on both the details of the SM-SC 
coupling at the interface and the dimensions of the nanosys- 
tem. More specifically, considering a SM wire with rect- 
angular cross section Ly x L^, the inter-band effects be- 
come critical whenever the gaps AE^ri' separating differ- 



ent confinement-induced bands n 



are compara- 



ble with the bulk SC gap Aq. This situation occurs in wide 
wires and in the "sweet spot" regime l^j^l, when two spin 
sub-band become degenerate at finite Zeeman field. By con- 
trast, in the limit AEnn' ^ Aq the inter-band effects are 
negligible and one can treat the proximity effect in the de- 
coupled band approximation, i.e., jnn' ~ if n ^ r^^ In 
this approximation, the matrices describing the proximity ef- 
fect become diagonal, e.g., jrin' = "yn^nn', and we have 
Zn l/Vl + 7n/Ao and A^ 7^Ao/(7r^ + Aq). Also, 
we note that there is a key difference between the proximity- 
induced inter-band effects associated with confinement in the 
directions normal and parallel to the interface, respectively. 
Specifically, assuming = n'^, i.e., strong confinement in 
the normal direction (so that the low-energy physics is con- 
trolled by a single mode), and a weaker confinement in 
the y-direction, along with with non-uniform SM-SC hop- 
ping i{iy), results in an effective SM-SC coupling jnn' = 
"jnyn' that contains non vanishing off-diagonal terms. How- 
ever, for typical SM-SC couplings, 7n^n' becomes negligi- 



ble when the difference In,, 



is large, as one can infer 



Oy ,Vy\ 

from Eq. ( [26| ), and only neighboring Uy bands are signifi- 
cantly coupled. Consequently, when constructing the effec- 
tive Hamiltonian ( [3T] ) it is enough to consider a reduced num- 
ber of Uy bands to obtain an accurate description of the low- 
energy physics. This scenario was investigated in detail in 
Ref. |70|. On the other hand, when the confinement in the 
z-direction is weak, bands with arbitrary become coupled, 
as 7n^n; ^ 4^nA^0z)<l^n'^{'^0z)' this casc, it is critical to 
include the high-energy bands in the calculation, as they 
re-normalize the low-energy physics via proximity-induced 
virtual processes. 



D. Low-energy physics in Major ana wires with disorder 

There are several types of disorder that may play a role in 
the low-energy physics of semiconductor-supraconductor hy- 



brid structures fTOl: impurities inside the bulk s-wave SC, dis- 
order in the SM nanowire, and disorder induced by random 
SM-SC coupling. Below, we mention a few representative 
sources that could generate these types of disorder and briefly 
discuss some of their main features. As a general remark, we 
point out that the theoretical treatment of disorder in semi- 
conductor Majorana wires can be done along two different 
directions: calculations of disorder-averaged quantities and 
calculations involving specific disorder realizations. How- 
ever, the low-energy physics of the nanowire (on a scale of 
the order of the induced gap) is controlled by a small number 
of quantum states, typically less than 100 (see, for example. 
Figs. [T6]and[T7] upper panels, in Appendix IdI) having most 



of their spectral weight inside the SM wire. When disorder 
is located in the nanowire or at the interface, the energies and 
the wave functions characterizing these states depend signifi- 
cantly on the specific disorder realization and the relevance of 
the disorder-averaged quantities is questionable. The reason 
for this dependence can be understood qualitatively by notic- 
ing that features of the disorder potential with length scales 
smaller than l/kp, where kp is the Fermi wave vector of the 
Majorana band, are irrelevant, since they are averaged over 
by the low-energy SM states. Typically, is of the order 
10^ nm, hence the effective disorder potential in a nanowire 
of length Lx ~ 1/im is characterized by a small number of 
scattering centers and, consequently, the details of the disor- 
der potential become relevant. Disorder-averaged quantities 
provide a good description of the low-energy physics (for ex- 
ample, the value of the induced SC gap) in long wires fTTS], 
but do not capture the specific properties of a given small seg- 
ment of that wire, or those of a short wire of similar length. 

Scattering off impurities inside a disordered bulk supercon- 
ductor has a negligible effect on the topological SC phase of 
the nanowire 1116111171 . In essence, this behavior is due the 
fact that the SM effective impurity scattering rate involves 
higher-order SM-SC tunneling processes and is suppressed 
by the destructive quantum interference of quasi-particle and 
quasi-hole trajectories 1 116|. Consequently, static disorder in 
the superconductor does not suppress the proximity induced 
topological superconductivity in the semiconductor. 

By contrast, disorder in the SM nanowire or at the SM- 
SC interface can strongly affect the stability of the topolog- 
ical SC phase. Some of the generic features of the low- 
energy states in the presence of disorder are illustrated in Ap- 
pendix |D] Possible sources of disorder are the random vari- 
ations of the width of the SM wire and, in general, surface 
roughness, and random potentials created by charged impu- 
rities located on or near the surface of the wire iTOl . An- 
other possible source of disorder is represented by the random 
coupling at the semiconductor-superconductor interface I70J . 
We note that this random coupling has a twofold manifesta- 
tion: a fluctuating induced pairing potential and a random 
proximity-induced renormalization of the SM Hamiltonian 
(see Sec. |IIIC|). In general, the suppression of topological 



superconductivity by disorder represents a serious challenge 
to the experimental realization of MFs. A recent proposal for 
optimizing the stability of the topological phase against dis- 
order involves replacing the semiconductor wire by a chain 
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of quantum dots connected by s-wave superconductors 11181 . 
In addition to the adverse effect on the stabiHty of the topo- 
logical SC phase, the presence of disorder impacts the low- 
energy physics of the SM nanowire-SC hybrid structures in 
a number of other ways, some of them being investigated 
theoretically in several recent studies. For example, it was 
shown that the topological quantum phase transition is char- 
acterized by a quantized thermal conductance and electrical 
shot noise power that are independent of the degree of disor- 
der |119|. The interplay between disorder and interaction in 
one-dimensional topological superconductors was addressed 
in Ref. 1 1201 . Also, it was shown that in systems with disor- 
der located at the end of the wire the weight of the characteris- 
tic Majorana-induced ZBCP in the differential tunneling con- 
ductance is strongly enhanced by mixing of sub-bands 1 121J. 
On the other hand, the presence of disorder can generate a 
ZBCP even when the superconducting wire is topologically 
trivial due to the proliferation of disorder-induced low-energy 
states 1 93 1, or due to the weak antilocalization resulting from 
random quantum interference by disorder 11221 . A similar 
ZBCP can occur in the topologically trivial phase as a result 
of fermionic end states with exponentially small energy that 
emerge if the confinement potential at the end of the wire is 
smooth | 92 |. In the light of the recent experiments on semi- 
conductor nano wire-superconductor hybrid structures, ruling 
out the alternative mechanisms for the conductance peak rep- 
resents a serious challenge. 



IV. RECENT EXPERIMENTS AND THEORETICAL 
INTERPRETATIONS 

In this section we briefly summarize several recent obser- 
vations of experimental signatures consistent with the realiza- 
tion of Majorana bound states in semiconductor nano wire- 
superconductor structures. We also discuss a number of ap- 
parent discrepancies between the observed features and the 
theoretical predictions based on simple model calculations for 
the Majorana wire. 



A. Experimental signatures of Majorana fermions in hybrid 
superconductor-semiconductor nanowire devices 

The observation of a zero bias conductance peak in lo- 
cal tunneling conductance measurements on semiconductor 
nanowires coupled to an s-wave superconductor has been re- 
cently reported in Ref. 1561 . This observation, which is con- 
sistent with the theoretical predictions, may represent the first 
experimental evidence of Majorana fermions in a condensed 
matter system. Soon after, observations of similar ZBCPs 
have been reported by two other groups fST, '581 • Since there 
are significant differences among the experimental setups, es- 
tablishing conclusively that the ZBCPs observed in differ- 
ent systems are due to the same mechanism remains a crit- 
ical open question. Below, we briefly summarize the main 
results of Ref. (561 . We mention that a measurement of 
the fractional a.c. Josephson effect has also been reported in 
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FIG. 7. (Color online) Top panel: Schematic representation of the 
experimental setup showing a cross-section through the device mea- 
sured in Ref. |56|. Bottom panel: Differential conductance dl/dV 
versus bias potential F at a temperature T = 70mK and for different 
values of the magnetic field from to 490 mT in of 10 mT. Traces 
are offset for clarity. Figure adapted from Ref. 1561 . 

Ref. 1591 . The observation suggests the presence of a Shapiro 
step with a height twice larger than the value expected for 
conventional superconductor junctions |59|, which would be 
consistent with the presence of MPs. 

As discussed in the previous sections, an optimal hybrid 
system capable of hosting Majorana bound states and of mea- 
suring the expected ZBCP should posses certain characteris- 
tics, such as i) a large SC gap, to protect the MF and to allow 
tuning the Zeeman field in a wide range without destroying 
the superconducting phase, ii) a long-enough wire length, to 
accommodate Majorana bound states at the wire ends without 
large overlap, iii) low disorder, iv) the ability to control the 
chemical potential, and v) the ability to tunnel from a normal 
lead, while minimizing the perturbation on the topological SC 
state. The basic experimental setup for the measurements re- 
ported in Ref. [56 1 consists of a InSb nanowire in contact with 
a superconductor (NbTiN) and a metallic (Au) lead, as shown 
schematically in the top panel of Fig. [7] To measure the differ- 
ential conductance dl/dV, a bias voltage is applied between 
the normal lead and the superconductor (SC) and a tunnel bar- 
rier is created in the region between the metallic lead and the 
SC segment of the nanowire by applying a negative voltage to 
a narrow gate. The InSb nanowire is characterized by a high 
^-factor (g ^ 50), as well as a strong spin-orbit (SO) cou- 
pling (the Rashba parameter is a ^ 0.2eV-A), making it a 
good candidate for the realization of MFs in the presence of 
proximity-induced superconductivity. The barrier potential 
suppresses the background conductance due to Andreev pre- 
cesses and helps to reveal the ZBCP generated by the possible 
presence of MFs that could emerge for values of the magnetic 
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FIG. 8. (Color online) Dependence of dl/dV on the magnetic field 
orientation. (A) Rotation in the plane of the substrate for a magnetic 
field |B| = 200 mT. Note that the ZBCP vanishes when B is perpen- 
dicular to the wire, i.e., parallel to the spin-orbit field. (B) Rotation 
in the plane perpendicular to the spin-orbit field for |B| = 150 mT. 
Note that a ZBCP is present for all angles. The top panels show 
line-cuts at angles marked with corresponding colors in (A) and (B). 
Figure adapted from Ref. l56l . 

field above a certain critical value. The dependence of the dif- 
ferential tunneling conductance on the bias voltage for differ- 
ent values of the magnetic field is shown in the bottom panel 
of Fig. |7] Note that a ZBCP emerges when the applied mag- 
netic field exceeds about 100 mT. The ZBCP is robust to fur- 
ther increasing the magnetic field until about 400 mT, above 
which the single peak seems to split into a two-peak structure. 
These observations are consistent with the MF interpretation 
of the tunneling conductance, as shown, for example, by the 
theoretical predictions illustrated in Fig. [4] 

An important consistency check on the MF interpretation of 
the ZBCP is the requirement that the peak must disappear if 
the angle between the applied magnetic field and the direction 
of the effective spin-orbit field vanishes. The likely orienta- 
tion of the spin-orbit field is along the direction perpendicular 
to the wire in a plane that is roughly parallel to the nanowire- 
SC and nano wire-substrate interfaces. The angle dependence 
of the ZBCP was measured in Ref. |56 | (see Fig. [8]) and found 
to be consistent with the Majorana scenario, thus strengthen- 
ing the identification of the observed ZBCP with the presence 
of a MF bound state at the normal lead-nanowire interface. 

The main observations reported in Ref |56J, in particu- 
lar the emergence of a ZBCP at finite magnetic field, have 
been subsequently confirmed by other groups |[57][58l. The 
similarities are rather surprising, considering that many of 
the relevant parameters that characterize the devices mea- 
sured in these experiments are very different from those of 
Ref [[56| . For example, the experiments reported in Ref. |58| 
use InAs nanowires in proximity to superconducting Al (in- 
stead of InSb and NbTiN). More importantly, the segment of 



the nanowire that is proximity coupled to the superconductor 
is very short (approximately 150 — 200nm, which is compara- 
ble to the SC coherence length). Is it possible to realize Majo- 
rana bound states in such a short nanowire? A possible answer 
is provided in Ref. 123 , which shows that, in a short wire, the 
lowest energy state represents a pseudo-Majorana mode that 
evolves continuously into a true zero-energy Majorana mode 
as the wire length increases. 

In the recently reported experiments, there are several strik- 
ing features, the most prominent being the absence of a sig- 
nature associated with the closing of the quasiparticle gap at 
the topological quantum phase transition (TQPT) that sepa- 
rates the topologically trivial and nontrivial SC phases. The- 
oretically, this vanishing of the gap at the critical magnetic 
field must precede the emergence of the Majorana-induced 
ZBCP ITT] [41] [TOl. This can be clearly seen in the simula- 
tion shown in Fig. |4| but is absent in the experimentally mea- 
sured tunneling conductance (see Fig. [8]). Another striking 
feature is the "soft" nature of the induced SC gap, as con- 
sistently revealed by the differential conductance measured in 
the recent experiments [56.-581. Furthermore, the observed 
ZBCP is more than an order of magnitude weaker than the 
quantized value of 2e^ /h predicted theoretically. Subsequent 
theoretical work f94l has attributed this discrepancy to a fi- 
nite temperature effect in conjunction with a hybridization- 
induced splitting of the Majorana mode in finite wires. In ad- 
dition to these unexpected experimental features, a number 
of alternative scenarios for the emergence of ZBCPs, which 
do not involve the presence of MFs but more conventional 
mechanisms involving strong disorder [i93lll22lll24l . smooth 
confinement |92|, or Kondo physics 11251 , have been recently 
proposed theoretically, making this problem an exciting but, 
at the same time, a rather fluid and confusing subject. Below, 
we discuss several theoretical proposals that provide possi- 
ble explanations for the experimentally observed features and 
suggest further tests that could better reveal the nature of the 
ZBCPs observed in charge transport measurements on semi- 
conductor nanowire-superconductor structures. We empha- 
size that the recent developments in this field reveal that key 
observable features are strongly dependent on various details 
of the system. Significant understanding of the relevant phys- 
ical mechanisms responsible for these features, at least at a 
qualitative level, can be achieved by realistically modeling the 
semiconductor-superconductor structures . 



B. Suppression of the gap-closing signature at the topological 
quantum phase transition 

While the observations (56I of a ZBCP in conductance 
measurements on semiconductor nanowires coupled to super- 
conductors may represent the first experimental evidence of 
MFs, the absence of any signature associated with the closing 
of the superconducting gap (see Fig. [7| at the critical mag- 
netic field asscociated with the topological quantum phase 
transition (TQPT) casts serious doubt on the MF interpreta- 
tion of the ZBCP. The closing of the quasiparticle gap at the 
TQPT is a fundamental theoretical requirement. The question 
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is whether or not this gap closure has a visible signature when 
the system is probed experimentally. Ref. 1 126] has offered 
a possible explanation for the observed non-closure of the 
gap at the quantum phase transition between the trivial super- 
conductor and the topological superconductor supporting the 
MFs in the end-of-wire tunneling experiments. By solving 
numerically an effective tight-binding model for multiband 
nano wires with realistic parameters, it has been shown |126| 
that, in the vicinity of the topological transition, the amplitude 
of the low-energy states at positions near the ends of the wire 
may be orders of magnitude smaller than the amplitudes of the 
localized MFs. Consequently, if the chemical potential of the 
system is located near the bottom of a confinement-induced 
semiconductor band, the contributions of these low-energy 
states to the end-of-wire local density of states (LDOS), and 
hence to the tunneling conductance, are essentially invisible 
(see Fig. |9] middle panel). This behavior results in an appar- 
ent non-closure of the gap, as reflected by these local quan- 
tities, even though the magnetic field approaches the critical 
value and the system is driven through a TQPT. By contrast, 
the closing of the gap mandated by the topological transition 
is clearly revealed by other quantities, such as the total den- 
sity of states (DOS) (see Fig. [9| top panel) and the LDOS 
near the middle of the wire (Fig. |9] bottom panel). A defi- 
nite prediction of Ref. 11261 is that a tunneling measurement 
near the middle of the wire should clearly reveal the closing 
of the gap at the critical field Fc, but the zero-bias peak asso- 
ciated with the presence of Majorana bound states should be 
absent. Correlated with an end-of-wire measurement charac- 
terized by a zero-bias peak for F > Fc, this would constitute a 
powerful argument for the presence of zero-energy Majorana 
bound states. In addition, the model calculations have shown 
that the non-closure of the gap, as revealed by the end-of- 
wire LDOS, is a non-universal phenomenon. For example, if 
a local measurement is performed in a regime characterized 
by a value of chemical potential that is not close to the bot- 
tom of a semiconductor band (i.e., A/i A), the closing of 
the gap should be visible. This regime is illustrated in Fig. [4] 
Finally, we note that the dominant feature in the end-of-wire 
LDOS (see Fig. |9] middle panel) are due to states associated 
with low-energy occupied bands that have significant ampli- 
tudes near the ends of the wire. In certain conditions (e.g., 
in the presence of disorder or in systems with smooth con- 
finement) these states can have energies lower than the bulk 
quasiparticle gap and, consequently, can provide substantial 
contributions to the in-gap LDOS. 



C. Discriminating between Majorana bound states and garden 
variety low-energy states 

Despite their conceptual simplicity, the zero bias conduc- 
tance peak experiments do not constitute a sufficient proof 
for the existence of Majorana bound states in semiconductor- 
superconductor hybrid structures. In a recent theoretical 
work (921, it has been shown that a non-quantized nearly 
zero bias peak, such as that observed in the recent experi- 
ments |[56ti58l . can arise even without end state MFs, pro- 
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FIG. 9. (Color online) Top: Density of states (DOS) as a function 
of B for a wire with four partially occupied bands (seven spin sub- 
bands) and a chemical potential near the bottom of the fourth band 
(for B = 0). At the TQPT characterized by Be ^ 0.2T the bulk gap 
closes. Middle: Local density of states (LDOS) at the end of the wire 
as a function of magnetic field. The strong finite energy features have 
a weak dependence on B. For B > Be the Majorana peak is clealy 
visible. Note that there is no visible signature of the bulk gap closing 
at the TQPT. Bottom: LDOS at the middle of the wire. Note the 
closure of the gap at the TQPT and the absence of the zero-energy 
Majorana peak. Figure adapted from Ref. LI 26 J . 



vided the confinement potential at the wire end is smooth. A 
non-quantized nearly zero energy peak at the wire ends has 
also been shown to occur due to strong disorder effects 1931 , 
even when the nano wire is in the topologically trivial phase. 
In essence, these nearly-zero bias peaks occur at finite val- 
ues of the Zeeman field due to the proliferation of low-energy 



states in the presence of disorder (see Sec. |III D| ) or in wires 
with smooth confinement (seeSec. |IVD| ) . To discriminate be- 
tween different possible mechanisms responsible for the zero 
bias conductance occurring in a tunneling experiment, a di- 
agnostic signature for the Majorana-induced ZBCP has been 
recently proposed |127 |. It was shown that, for smooth con- 
finement at the ends of the wire, the emergence of the near 
ZBCPs is necessarily accompanied by a signature similar to 
the closing of a gap in the end-of-wire local density of state. 
This signature occurs even though there is no corresponding 
quantum phase transition, as the system stays in the topolog- 
ically trivial phase, and traces the Zeeman field dependence 
of the nearly-zero energy states. In the absence of such a gap 
closing signature, a ZBCP is unlikely to result from the soft 
confinement effect L92J . Similarly, when the ZBCP appears 
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at the wire ends from disorder effects (but without MFs), the 
emergence of the zero conductance is preceded by a signa- 
ture similar to the closing of the gap 193111271 . So far, among 
all the scenarios that have been considered, the topological 
phase transition scenario involving the emergence of the MFs 
is the only one consistent with a ZBCP that occurs beyond a 
certain critical magnetic field and the apparent non-closure 
of the quasiparticle gap before the emergence of the ZBCP. 
Since this is precisely what is observed in the experiments in 
Ref. f56|, these theoretical results strengthen the identification 
of the observed ZBCP with topological Majorana bound states 
localized at the ends of the wire. 

We note that the presence of low-energy sub-gap states in 
systems with disorder and smooth confinement may have an- 
other important consequence with dramatic experimentally- 
observable implications: the existence of a soft superconduct- 
ing gap, i.e., a SC gap characterized by a non-vanishing den- 
sity of states and a v-shaped sub-gap tunneling conductance 
even at very low temperatures. This feature is manifestly 
present in the experimental data 1561 (see, for example. Fig. 
??) and represents one of the most important open issues in 
this field. There are several critical aspects that need to be 
clarified, such as the origin of the in-gap spectral weight, why 
these in-gap contributions occur as a smooth background, 
rather than sharply defined peaks in the LDOS, and whether 
or not a Majorana-induced ZBCP can be well-defined in a 
system with soft gap. Recently, disorder induced by interface 
fluctuations was identified 11281 as the likely source of the 
in-gap states responsible for the soft gap. Another theoreti- 
cal work 1 127] has shown that there are two key ingredients 
that may explain the emergence of the soft gap in weakly con- 
fined wires: i) A finite potential barrier allows states with large 
spectral weight near the end to hybridize with metallic states 
from the leads, ii) States from lower-energy bands can pene- 
trate through the barrier, hybridize strongly with the metallic 
states, and generate broad contribution to the LDOS. By con- 
trast, the Majorana mode, which is associated with the top 
occupied band, couples weakly to the lead, hence it is weakly 
broadened and still generates a well defined ZBCP on top of 
the smooth background. 



Majorana physics in finite-size wires: A "smoking gun" for 
the existence of the Majorana mode 



D. 



It has recently been proposed flOO] that direct observation 
of the splitting of the zero bias conductance peak could serve 
as a "smoking gun" evidence for the existence of the Majo- 
rana mode. In essence, the Majorana bound states come al- 
ways in pairs | 23 1 that are localized near the ends of the wire, 
if the system is clean enough. In any finite wire, the wave 
functions of the two Majoranas overlap, leading to a split- 
ting 1 94, 123 ,fT29l[T30| of the Majorana mode, which is a 
pure zero-energy mode only in the infinite wire limit. At 
fixed chemical potential, this hybridization-induced energy 
splitting is characterized by an oscillatory behavior that de- 
pends on the Fermi wave vector of the top occupied band and 
the length of the wire. By contrast, when the particle den- 
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FIG. 10. (Color online) Semiconductor nanowire with four partially 
occupied bands in the presence of a nonuniform potential V{x). The 
position-dependence of V is shown in panel (a). The profiles of sev- 
eral low-energy states with \Ei\ < \E2\ < ... correspond- 
ing to different values of the Zeeman field are shown in panels (b) 
r = 0.25meV (trivial SC phase) , (c) F = 0.4meV (topological SC 
phase), and (d) F = O.SmeV. Note that the wire is effectively sep- 
arated into two segments for states with a small characteristic wave 
vector. States corresponding to large kp values penetrate through 
the potential barrier AV. With increasing Zeeman field, low-energy 
states associater with the low-energy bands emerge because of the 
smooth confinement, e.g., the state n = 3 (green line) in panel (d). 



sity is constant, the oscillations can be suppressed flOO], but 
the splitting of the Majorana mode is still a generic feature. 
Furthermore, regardless of conditions, two independent tun- 
neling measurements at the opposite ends of a wire should 
observe exactly the same splitting of the ZBCP [ 1QQ |, as long 
as: i) the peak is due to Majorana splitting, and ii) a single 
pair of Majorana bound states exists in the system. We em- 
phasize that features other than the splitting itself observed 
in local measurements at the two ends of the wire may be 
completely different. If confirmed experimentally, this would 
constitute strong evidence for the existence of the elusive Ma- 
jorana mode in semiconductor-superconductor structures. 

In the presence of disorder, several states with energies 



lower than the induced SC gap will emerge (see Sec. HID). 
The energies of these sub-gap states depend on the strength 
of the disorder and the values of the Zeeman field. However, 
the Majorana mode is still protected, as long as the disorder 
strength does not exceed a certain critical value, and, conse- 
quently, the hallmark signature described above should still 
be observable. Nonetheless, strong disorder can effectively 
cut the nanowire into two or more segments that can each 
host a pair of Majorana bound states. A natural question con- 
cerns the fate of the splitting oscillations in the presence of 
a strong scattering center that could effectively cut the wire 
into two segments, each hosting a pair of Majorana bound 
states. To address this questions, we consider we consider 
a wire with rectangular cross section Ly x = 100 x 40 
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FIG. 11. (Color online) Same as in Fig. [To] but for a lower value of 
the potential barrier AV. Note that, unlike Fig. [To] the lowest energy 
state (n = 1, red line), corresponding to two overlapping Majorana 
bound states, is peaked near the ends of the wire for Zeeman fields 
associated with the topological SC phase. 

nm, four palatially occupied bands, and smooth confinement 
in the presence of a strong scattering center consisting in a po- 
tential barrier of hight AV located at distances L^i ^ 0.6/im 
and Lx2 ^ 1.2/im from the two ends of the wire, respec- 
tively. The confinement and the barrier are provided by the 
position-dependent potential V{x) shown in the upper pan- 
els of Fig. [To] and Fig. [TT] When F > = 0.35 meV 
and AV = 0, two Majorana bound states are localized near 
the ends of the wire (i.e., x ^ 0.2/im and x ^ 2.2/im, re- 
spectively), while for large AV^ the wire is split into two dis- 
connected segments, each of them hosting a pair of Majo- 
ranas. The behavior of the low-energy states in the interme- 

As a general 



diate regime is illustrated in figures and 1 1 



feature, we note that states corresponding to small character- 
istic wave vectors (i.e., long wavelength oscillations) tend to 
be contained in one of the two segments of the wire separated 
by A V. By contrast, states with large characteristic wave vec- 
tors (rapid oscillations) can more easily penetrate through the 
finite barrier. Furthermore, we note that Fermi /c-vector kp 
associated with the top occupied band increases with the Zee- 
man field, while the Fermi wave vectors corresponding to the 
lower energy bands are always larger than kp- Consequently, 
in the presence of a large barrier, the two pairs of Majoranas 
are characterized by a strong intra-pair hybridization and a 
weak inter-pair hybridization, which leads to two low-energy 
modes mostly localized inside the L^i and segments, re- 
spectively. The corresponding wave functions (states n = 1 - 
red line - and n = 2 - yellow) are shown in Fig. [T0| Reducing 
AV (or increasing the Zeeman field) increases the inter-pair 
hybridization, which results in a lowest energy state with max- 



ima near the ends of the wire (see Fig. 11). 

The hybridization-induced splitting oscillations of the Ma- 
jorana mode can be clearly seen in the field dependence of 
the density of states (DOS), as shown in Fig. [12] We 
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FIG. 12. (Color online) Density of states as a function of the Zeeman 
field in the presence of a nonuniform potential with AV — ImeV 
(top) and AV = 2meV (bottom). Note the closing of the bulk gap at 
Fc ~ 0.35meV and the oscillations of the lowest energy modes for 
r > Fc. The corresponding LDOS at the opposite ends of the wire 
is shown in Fig. [13] 

want to emphasize two features. First, regardless of the 
number of coupled Majorana pairs, only one mode exhibits 
zero-energy crossings at discrete values of the Zeeman field, 
which is a characteristic signature of Majorana physics in 
finite nanowires il23l . By contrast, the other low-energy 
mode is characterized by minima that vanish only in the limit 
AV oo, i.e., vanishing inter-pair tunneling. Second, the 
period and the amplitude of the oscillations induce by intra- 
pair hybridization increase with decreasing wire length ap- 
proximately as l/Lx IIOOI . Since Lx2 ^ 2La^i, the ratio 
between the corresponding periods and amplitudes is approx- 
imately two, as evident from the lower panel of Fig. [12] We 
note that, in very short wires (quantum dots), the period of the 
oscillations may be large, so that practically only one zero- 
energy crossing is accessible |123|. Also, in such quantum 
dot-superconductor structures the Coulomb interaction is ex- 
pected to play an important role 1 131 1. 

The characteristic signatures of the splitting oscillations in 
a set of two independent local measurements at the opposite 
ends of the wire, as reflected by the corresponding local den- 



sity of states (LDOS), are shown in Fig. 13 Several features 



need to be emphasized. As discussed in Sec. IV B the sup- 
pression of the gap closing signature stems from the spatial 
properties of various low-energy states. In short wires, there 
is no qualitative difference between extended states and states 
localized near the wire ends, and consequently, a gap closing 
signature should be visible. This signature should be stronger 
at the left end of the wire, since the segment Lxi is shorter. 
However, the corresponding features are rather weak (see Fig. 
13 ) and could be very hard to resolve when on top of an in- 
coherent in-gap background, such as that responsible for the 
soft gap observed experimentally |56-58|. In an actual exper- 
iment, the low-energy states couple to metallic states from the 
leads, while the states associated with low-energy bands are 
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FIG. 13. (Color online) Splitting oscillations of the ZBCP as re- 
flected by the LDOS at the left (L) and right (R) ends of the wire. 
The top two panels correspond to the parameters of Fi g. pT] while 
the lower two panel are for the system illustrated in Fig. |10[ 

expected to have significant weight outside the superconduct- 
ing region of the wire \127] . This leads to strong broadening 
and could explain the the soft gap observed in the experiments 



while hiding all the weak features that are present in Fig. [T3_ 
Consequently, we expect the splitting oscillations measured at 
the two ends of a wire containing a strong scattering center to 
be uncorrected, as shown by the dominant low-energy fea- 



tures corresponding to Ay = 2meV in Fig. [13] By contrast, 
in the presence of a weak scattering center (AV = ImeV in 
Fig. [13]), the dominant low-energy features reveal the same 
ZBCP slitting at both ends of the wire. However, the ampli- 
tude of the splitting is strongly reduced because of the larger 
Lx and the oscillations may not be visible at finite energy res- 
olution. Therefore, we conclude that the optimal experimental 
setup for observing correlated Majorana-induced splitting os- 
cillations should use high quality wires with < Ifim. 



V. CONCLUSIONS 

The emergence of a zero-bias conductance peak above 
a certain critical value of the Zeeman field represents the 



necessary condition for the existence of zero-energy Ma- 
jorana bound states in spin-orbit-coupled semiconductor- 
superconductor hybrid structures. Strong experimental evi- 
dence for signatures consistent with Majorana physics were 
recently reported in charge transport measurements on such 
nanowire heterostructures. These encouraging experimental 
developments have also raised a number of rather unexpected 
questions that need to be clarified before claiming victory. 
Nonetheless, corroborating these observations by performing 
a "smoking gun" measurement, such as the proposed obser- 
vation of the correlated splitting of the zero-bias peak, could 
be a step towards unambiguously establishing the existence of 
Majorana bound states in semiconductor nano wires. This goal 
could be achieved in the near future. In addition, a fractional 
AC Josephson effect, 2e^ /h quantized conductance through 
a MF state, or a signature of a Zeeman-tuned TQPT, say, by 
tunneling to a region away from the wire ends, can decisively 
establish the existence of MFs in semiconductor nanowire het- 
erostructures. 

However, the search for the non-Abelian Majorana fermion 
zero modes in semiconductor heterostructures is far from over. 
Future experiments have to address the sufficient conditions 
for the existence of the Majorana mode, which may consist 
a tunneling based interference measurement such as that dis- 
cussed in Sec. IIH or the direct observation of non-Abelian 
braiding statistics in some form. This direction will involve 
challenging experimental problems regarding the controlled 
engineering of complex devices, as well as basic aspects re- 
lated to quantum decoherence and the manipulation of Majo- 
rana bound states. An important lesson provided by the re- 
cent developments is that, in real systems, key observable fea- 
tures are determined by various details of the structure, despite 
the topological nature of the superconducting state predicted 
to host the Majorana bound states. Therefore, to be able to 
discriminate between alternative scenarios and to clarify the 
mechanisms responsible for various features observed in the 
experiments, it is critical to develop realistic models of the 
heterostructures as discussed in Sec. III. 
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Appendix A: Inter-band pairing and the "sweet spot" 

In the 'standard model' of the Majorana nanowire, which 
involves a semiconductor (SM) wire with spin-orbit coupling, 
proximity-induced superconductivity, and Zeeman spin split- 
ting, the emergence of zero energy Majorana bound states re- 
quires values of the chemical potential consistent with an odd 
number of partially occupied sub-bands. As an example, let 



us consider the case illustrated in Fig. 14 involving two bands, 
n = 1 and n = 2, and four values of the Zeeman field F. The 
splitting between the pairs of sub-bands with opposite helicity 
(i.e., n— and n+) is proportional to the applied field and, for 
r = Fq, the sub-bands 1+ and 2— are degenerate at /c^^ = 0, 
ei+(0) = 62- (0) = /i2. The condition for the existence of 
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FIG. 14. (Color online) Schematic representation of the evolution 
of the SM nanowire spectrum with the applied Zeeman field. The 
curves represent the energy as function of the wave vector kx for 
four different values of the Zeeman field F. For F = Fo (not shown) 
the sub-bands 1+ and 2— are degenerate at /c^ = 0. In the presence 
of induced superconductivity, a system with a value of the chemical 
potential /x = in or fi — /xs is in a topological SC state, while 
for 11 = the SC state is topologically trivial. However, in the 
presence of inter-band paring, A 12 / 0, the SC state corresponding 

10 11 — 112 becomes topologically nontrivial in the vicinity of F = Fo 
(the sweet spot). 

a topological superconducting (SC) phase is satisfied for all 
four values of F when the chemical potential is /i = /ii (one 
partially occupied sub-band) or /i = /ia (three partially occu- 
pied sub-bands), but the system is topologically trivial when 

11 = 112 (two occupied sub-bands). This picture holds as long 
as the inter-band paring is zero. However, when A12 ^ 0, 
the number of partially occupied sub-bands does not repre- 
sent a good criterion for establishing the topological nature 
of a given SC phase. More specifically, a quasi- ID nanowire 
with inter-band pairing and control parameters in the vicinity 
of the "sweet spot" F — Fq and /i = /i2 is in the topologi- 
cally nontrivial SC phase provided the inter sub-band spacing 
1^1+ — 62- 1 is smaller that the off-diagonal coupling A12 1691 . 
In the limit of strong inter-band mixing, the sweet spot regime 
is characterized by a topological state that is robust against 
chemical potential fluctuations, such as those created by dis- 
order. This stronger immunity can be understood in terms of 
the range of chemical potentials consistent with a topological 
SC state for a given value of the Zeeman field. Considering 
example shown in Fig. 14 for F = 0.4Fo there are two narrow 
topological regions (for chemical potentials in the vicinity of 
/ii and /i3, respectively), separated by a topologically trivial. 
By contrast, for F sufficiently close to Fq the range of ji con- 
sistent with a topological SC state extends from the bottom of 
sub-band 1— to the bottom of 2+, without passing through 
any topological quantum phase transition. 

How can inter-band mixing occur in SM nanowires with 
proximity-induced superconductivity? In essence, non- 
uniform interface tunneling generically induces inter-band 



coupling. Consider, for example, a hybrid structure contain- 
ing a SC and a long SM nanowire with rectangular cross sec- 
tion with Ly Lz. The low-energy physics is controlled 
by the lowest energy transverse mode, Uz = I, but , in gen- 
eral, involves several Uy bands. Now let us assume that only 
half of the nanowire is covered by the superconductor, so that 
i{y) = ifor y < Ly/2 and i{y) = Ofor y > Ly/2. Using Eq. 
([26|), the effective SM-SC coupling becomes 



Ny/2 

Inyn'y =70 ^ Ky(iy)(t)n'y(iy), 



(Al) 



where 70 = ^f|^P [^1(^02)]^- The diagonal coupling is 
luyUy = 7o/2, but, in addition, each band couples to all 
other bands of opposite parity. The inter-band couplings for 
neighboring bands, n'y = Uy ± 1, are comparable to the di- 
agonal value, e.g., 712 ~ 0. 85711. As described in Section 



IIIC[ the off-diagonal SM-SC coupling generates proximity- 
induced inter-band pairing, A^^n' , which leads to the sweet 
spot physics described above. Finally, we note that the exper- 
imental realization of non-uniform SM-SC coupling by par- 
tially covering the SM nanowire with superconductor 1 56 1 can 
provide additional benefits, such as reducing SC-generated 
screening and allowing gate-voltage control of the chemical 
potential. 



Appendix B: Projection of the effective tight-binding 
Hamiltonian onto a low-energy subspace 

To realize the projection of the effective tight-binding 
Hamiltonian onto a low-energy subspace, one needs to iden- 
tify a convenient low-energy basis. While, in general, this 
problem has to be addressed using a combination of analyt- 
ical and numerical tools, the 2-band lattice model described 
above has a simple analytical solution. Specifically, for an 
electron-doped nanowire with rectangular cross section and 
dimensions ^ Ly Lz, the single particle quantum 
problem corresponding to from Eq. ([18]) has eigenstates 

V^na(*) = Ha^I ^nA(u)Xa, whcrC Tl = {fix, fly, Uz) with 

1 < < ^A. Xcr is an eigenstate of the spin operator, and 



sm ■ 



TrnA^A 



(Bl) 



In Eq. ( |B1| ) A = x^y^z and Lx = aNx, where a is the lattice 
constant. The eigenvalues corresponding to ipria are 



-2to cos 



TTTlx 



TTUy 



TTUz 



-cos ,^ ^ +COS — 6 I — /i, 

(B2) 

where ji is the chemical potential. 

Next, we assume that only a few bands are occupied, and 
that the low-energy subspace is defined by the eigenstates sat- 
isfying the condition Cn < Cmax, where the cutoff energy emax 
is typically of the order lOOmeV. Using this low-energy basis, 
the matrix elements of the total Hamiltonian can be written 
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explicitly, as described in Ref. fTOl. The matrix elements of 
the SOI Hamiltonian from Eq. ([18]) are 



Uzu'^ I Q.nxu'^^'^^y) (T(t' ^UyT 



where 



1 _ (^_l)nA+n';^ 



Nx + l 



COS 



(B3) 



(B4) 



cos 



Nx + l 



The first term in Eq. ( |B3| ) represents the intra-band Rashba 
spin-orbit interaction, while the second term couples differ- 
ent confinement-induced bands. Similarly, assuming that the 
Zeeman splitting V is generated by a magnetic field oriented 
along the wire (i.e., along the x-axis), T = g* jisBx/'^, where 
^* is the effective g-f actor for the SM nanowire, the matrix 



elements for the corresponding term in Eq. ( 15 ) are 



where a = —a. Adding together these contributions, the ef- 
fective Hamiltonian describing the low-energy physics of the 
semiconductor nanowire in the presence of a Zeeman field be- 
comes Hnn^ = {tpn\HsM + i^z|V^^), with tpn representing 
the spinor {ipn^j i^ni)- Explicitly, we have 



(B6) 



where n = (n^^, n^, n^), Cn is given by Eq. (B2) and Quxn'^ 
by Eq. ( |B4[ ). Note that a similar effective low-energy Hamil- 
tonian can be written for an infinite quasi- ID wire. In the 
limit Lx ^ oo, the wave vector kx becomes a good quantum 
number and we have 



Hnn'{kx) = [en{kx)^OiRkxCry^Tax]Sr, 



'^^Q.nyn'y^x^nzn'^ •> 

(B7) 



where = aa and n = {ny^Uz) labels the confinement- 
induced bands with energy 



(^n{kx) = 1^ - 2^0 ( cos ■ ^ 



2m* 



-cos 



_nriz_ 
7V, + 1 
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FIG. 15. (Color online) Dependence of the effective SM-SC cou- 
pling of the top valence band in a hole-doped SM wire with rectan- 
gular cross section Ly x Lz on the Ly/Lz ratio. The strength of the 
coupling decreases with increasing the wire thickness Lz but, in con- 
trast with electron-doped wires, it shows a strong dependence on the 
wire width Ly. In the limit Ly ^ oo the coupling vanishes because 

(B5) 

the top valence band has purely heavy-hole character, i.e., = 0. 



of the localized states that define the effective SC Hamiltonian 
(22). Assuming that jj^sc lies within a band with s-type charac- 
ter, we notice that the orbitals that are responsible for the cou- 
pling across the interface are the \Z) orbitals. Consequently, 
only the eigenstates couple to the SC. The corresponding 
matrix elements in equation ([23]) can be written as 



*OjO 



2^ (^io+d jo^sai 



(CD 



where s\ 



and t is a constant, while the coupling of 

the 3/2 states vanishes, t- ^. =0. We address the follow- 
ing question: What is the consequence of this selective cou- 
pling on the strength of the proximity effect in hole-doped 
nanowires with rectangular cross section? Specifically, we fo- 
cus on the effective SM-SC coupling 7^ for the top valence 
band. 

To determine the effective coupling, we calculate numer- 
ically the wave functions ^± corresponding to the dou- 
ble degenerate top valence band at kx = 0. These states 
can be expressed as four-component spinors, '^si^^'^z) = 
[?/^ , V^i^l , V^i^s and we have 



Appendix C: Proximity effect in hole-doped semiconductors 

To illustrate this the dependence of the SC proximity effect 
on the details of the SM-SC coupling, we consider the case 
of hole-doped SM nanowires. We use the 4-band Luttinger- 



type model defined by equations (20) and (21 ) to describe a 
hole-doped SM nanowire with rectangular cross section prox- 
imity coupled an s-wave SC. We assume that the z axis is 
perpendicular to the interface and corresponds to the (0, 0, 1) 
crystal axis of the underlying fee lattice. The specific form of 
the SM-SC coupling Hamiltonian ([23]) depends on symmetry 



^F|t?^|V^^l(z„io.)|'. (C2) 



For comparison, in the case of an electron-doped wire with 
uniform tunneling across the SM-SC interface, the effective 
coupl ing of the lowest conduction band can be obtained from 
Eq. ([26) for i{iy) = i md n = n' = (1,1,1). We have 
7c = i^Flil^^liioz) - 2n^VF\i?/{N, + 1)^. Note that 7e 
is independent of the wire width Ly and decreases approxi- 
mately as 1/L1 with increasing wire thickness L^, due to the 
decrease of the wave function amplitude at the interface. By 



24 



> 


300 




150 













c 


-150 


W 


-300 



r=0.6meV 
A|Li=0 






o 


o 

■ o o 




•v,=o 

°Vd=5meV 



-20 -10 10 

State index (n) 



20 



150 
100 
50 




bo o OqO 
oo 

Or^ 



300000 



098 

r^O O^ 



oooo 



OOq b o» ^o^ 

o ^00,00^8800898^0^ 

o -°o- 

o o o o o©^ 



1 2 3 4 5 6 
Disorder strength (meV) 



> 

D- 



400 
200 


-200 
-400 





r=0.6meV 
A(j=3meV 




,oO»- 

00 


oo 

:..«oo» 




•v,=o 

°Vd=10meV 



-20 -10 10 

State index (n) 



20 



> 



200 
150 
100 
50 




008000000, 



?88°°°ooo°o 



8888§°oooooooo° 

oo^ 



"Oo. 00^ 



'Oo o o 

oooooooooooooooo^ Oo^o^o,^^ 

*^'^Oo'^Oo°'^^888oooooooooo 
Oq o^ oooooooooo 
^o^Oo 
o^ Oo 
o Ooo^ 
Ooo^ "^oooooo 



2 4 6 8 10 12 
Disorder strength (meV) 



FIG. 16. (Color online) Top: BdG spectrum of a SM wire with 
proximity-induced superconductivity and Zeeman splitting F = 
0.6meV. The chemical potential relative to the bottom of the third 
band is A/i = 0. In a clean wire (Vd = 0) the Majorana bound states 
(red diamonds) are protected by a gap of the order ISO/xeV. This gap 
collapses in the presence of strong disorder. Bottom: Dependence of 
the non-zero lowest-energy modes on the strength Vd of the disorder 
potential. 



FIG. 17. (Color online) Top: BdG spectrum of a SM wire with pa- 
rameters corresponding to the topologically trivial phase. Note that 
there are no zero-energy states. The quasiparticle gap collapses in 
the presence of strong disorder. Bottom: Dependence of the lowest- 
energy modes on the strength Vd of the disorder potential. Note that 
the proliferation of disorder-induced low-energy states requires a 
stronger disorder potential that in the topological SC phase (see Fig. 




Appendix D: Low-energy states in the presence of disorder 



contrast, the effective coupling of the top valence band is char- 
acterized by a strong dependence on Ly for a fixed wire thick- 



ness Lz. The numerical results are shown in Fig. 15 Note 



that in the quasi-2D limit, Ly ^ 00, the effective coupling 
of the top valence band to the bulk SC vanishes. This can be 
understood by noticing that the top valence band of a SM slab 
is a purely heavy-hole band and, for a (0, 0, 1) surface orien- 
tation, the corresponding states have vanishing iIj^^I compo- 

nents. Consequently, we conclude that no proximity-induced 
superconductivity can occur in the top valence band of planar 
SM-SC heterostructures with a (0,0,1) interface. This re- 
sult holds for different interface orientations. Nonetheless, in 
quasi-lD nanowires, the top valence band is a superposition 
of heavy-hole and light-hole states and, consequently, the ef- 



fective coupling to the SC is nonzero, as shown in Fig. 15 



Note that the amplitude of the of the light-hole component 
varies non-monotonically with Ly and has a sharp disconti- 
nuity for a width-to-thickness ration Ly/L^ ~ 2. A quan- 
titative description of the SC proximity effect in hole-doped 
nanowires, requires a more detailed modeling of the wire (e.g., 
using the 8-band model). 



We illustrate some of the generic features of low-energy 
BdG spectrum of a disordered SM-SC hybrid system by con- 
sidering a SM nanowire of length Ix = 3/im and rectangu- 



lar cross section with L,, 



80nm and Lz 



40nm in the 



presence of a disorder potential V(r) = Vdfv{r), where the 
disorder profile is described by the random function fv with 
< 1 and Vd represents the amplitude of the disorder 
potential. The low-energy spectrum is calculated using a 2- 
band model, as described in Sec. |III A| All the calculation 
are done for a specific disorder realization, i.e., a fixed profile 
fv, but for variable disorder strength Vd- First, we consider 
a system with three partially occupied bands (five spin sub- 
bands) in the presence of a Zeeman splitting F = 0.6meV and 
a fixed chemical potential A/i = 0, as measured relative to 
the bottom of the third band in the absence of Zeeman split- 
ting. The clean wire (Vd = 0) is in the topological SC phase 
that supports zero-energy Majorana bound states. The corre- 
sponding BdG spectrum is shown in the top panel of Fig. 16 
(black circles and red diamonds for the Majorana states). In 
the presence of disorder (Vd 7^ 0), the mini-gap that protects 
the Majorana bound state becomes smaller and, eventually, 
collapses. The dependence of the non-zero lowest-energy 
modes on the strength of the disorder potential is shown in 
the bottom panel. We note that, in the presence of disorder. 
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FIG. 18. (Color online) Position dependence of the amplitude of the 
lowest-energy states for a wire with Lx = 3 /am and F = 0.6meV. 
The top panel corresponds to a clean wire in the topological SC 
phase. The red line with maxima near the ends of the wire repre- 
sents the Majorana bound states. In the presence of disorder, all the 
states become localized in various regions of the wire. The wave 
function of the lowest energy state of clean wire (Vd = 0) in the 
topological trivial state is characterized by an envelope with max- 
ima near the ends of the wire and oscillations corresponding to a 
certain value of kp- The lowest-energy states associated with the 
low-energy occupied bands (not shown) have similar characteristics, 
but higher values ofkp- For given values of the disorder strength and 
Zeeman splitting, states with larger characteristic Fermi wave vectors 
(shorter oscillation period) exhibit weaker localization. 



multiple zero-energy modes are possible as disorder effec- 
tively cuts the wire into disconnected topologically nontriv- 
ial segments supporting Majorana bounds states at their ends. 
Typically, these segments are relatively short and the states lo- 
calized at the ends of each segment (which can be viewed as 
representing a Majorana chain) overlap significantly and are 
characterized by nearly-zero energies that oscillate with the 
chemical potential and the Zeeman splitting. 

A similar proliferation of low-energy states with increas- 
ing the strength of the disorder can be seen in a system with 
parameters corresponding to the trivial SC phase, as shown in 



Fig. [17] Note that in this case the decrease of the quasiparticle 
gap with Vd is slower than that corresponding to A/i = 0. This 
behavior can be understood qualitatively by noting that disor- 
der tends to localize the low-energy states and that this effect 
depends of the characteristic Fermi wave vector associated 
with those states. Increasing A/i corresponds to larger values 
of kp and to a weaker localization. This behavior is illustrated 



by the profiles of the low-energy states shown in Fig. 18 The 
low-energy states corresponding to A/i = exhibit stronger 
localization that those for A/i = 3meV, in spite of the smaller 
amplitude of the disorder potential. We note that the low- 
energy states associated with the low-energy bands are char- 
acterized by large values of kp and, consequently, are harder 
to localize. Finally, we note that the other key parameter that 
controls the effect of disorder on the low-energy states is the 
Zeeman splitting. Increasing F, which breaks time-reversal 
symmetry, facilitates the collapse of the SC gap. Again, this 
effect is stronger for states with low values of the characteris- 
tic /ci?, which typically corresponds to the top occupied band, 
although significant band mixing is possible for certain types 
of disorder. Consequently, in a wire with hard confinement 
the nearly-zero energy states with most of the spectral weight 
coming from low-energy bands require higher values of the 
Zeeman splitting than the disorder-induced nearly-zero en- 
ergy states associated with the top band. On the other hand, 
the energy of these states can become exponentially small 
even in the absence of disorder if the confinement potential 
is smooth 1921 . Furthermore, the states associated with low- 
energy bands can penetrate through finite barrier potentials 
and extend into the normal section of the wire LI 27 J . 
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